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The theory of perturbations of spherically symmetric self-similar black holes is presented, in the 
Newman- Penrose formalism. It is shown that the wave equations for gravitational, electromagnetic, 
and scalar waves are separable, though not decoupled. A generalization of the Teukolsky equation 
is given. Monopole and dipole modes are treated. The Newman-Penrose wave equations governing 
polar and axial spin-0 perturbations are explored. 
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I. INTRODUCTION 



It has been known for decades that something strange and dramatic must happen near the inner horizon of a 
realistic black hole. Penrose (1968) [ll] first pointed out that a person passing through the Cauchy horizon (the 
' outgoing inner horizon) of an empty charged (Reissner-Nordstrom) or rotating (Kerr-Newman) black hole will see 
• the outside universe infinitely blueshifted. Empty here means empty of charge, matter, or radiation (but not of a 
r— 1' static electric fieldjexcept at the singularity. Subsequent perturbation theory investigations, starting with Simpson 
O & Penrose (1973) [86] and culminating with Chandrasekhar & Hartle (1982) '2(?|, confirmed that an observer crossing 
^ ' the Cauchy horizon of the Reissner-Nordstrom geometry will experience an infinite flux of radiation from the outside 
5-H , world. The prediction that an infinite energy density will collect at the Cauchy horizon is at odds with the assumption 
, that the black hole is empty, except in the exceptional case of an uncharged, non-rotating (Schwarzschild) black hole, 

whose inner horizon coincides with the singularity at zero radius. 
■ The full nonlinear nature of the instability near the inner horizon was clarified in a seminal paper by Poisson & 
K>" ' Israel (1990) [S^. Poisson & Israel argued that, if ingoing and outgoing streams are simultaneously present just 
0^ I above the inner horizon of a spherical charged black hole, then relativistic counter-streaming between the ingoing and 
outgoing fluids will lead to an exponentially growing instability which they dubbed "mass inflation" . During mass 
, inflation, the interior, or Misner-Sharpe [59|, mass, a gauge-invariant scalar quantity, exponentiates to huge values. 
• ■ Other gauge-invariant measures, such as the proper density, the proper pressure, and the Weyl scalar, exponentiate 
' along with the interior mass. The phenomenon of mass inflation has been confirmed analytically and numerically in 
^ . many papers S i, i, [11 Hi Efl El III • 

I Real black holes probably rotate (and rotate rapidly), but probably have very little charge, thanks to the huge 
• • . charge-to-mass ratio of individual protons and electrons, e/m-p « 10^*, where e is the dimensionless charge of the 
proton or electron, the square root of the fine-structure constant, and m-p is the proton mass in units of the Planck 
mass. The literature on mass inflation in rotating black holes is exceedingly modest [zl, [t^ (see also the review [5]), 
\^ • and there is much room for further investigation of this hard problem. In the meantime the present paper follows the 
5^ ] bulk of the literature in assuming that charge is a satisfactory surrogate for angular momentum (e.g. [24|). Physically, 
the negative pressure, or tension, of the radial electric field of a charged black hole produces a gravitational repulsion 
analogous to, and perhaps an adequate model for, that produced by the centrifugal force in a rotating black hole. It 
is this gravitational repulsion that is responsible for the presence of inner horizons inside empty black holes, within 
which the fiow of space becomes subluminal. 

Physically, what causes mass inflation? It is clear that relativistic counter-streaming between ingoing and outgoing 
streams is a central ingredient, as originally suggested by Poisson & Israel [83|. Hansen et al. ^l| have suggested 
that adiabatic compression may be a contributing factor. In our own previous investigations we considered feeding a 
spherical black hole in a self-similar fashion both with a single fluid of charged baryons [39] , and with two counter- 
streaming fluids, charged baryons and neutral dark matter [40|. We found that mass inflation did not occur in the 
case of a single fluid, but did occur with two fluids. For the baryons we assumed a relativistic equation of state 
p — ip, so that the sound speed, at a/1/3 of the speed of light, was less than the speed of light. Thus in the case of 
the single fluid there was no possibility of relativistic counter-streaming, and indeed no mass inflation occurred. In 
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FIG. 1: Partial Penrose diagram of the Reissner-Nordstrom geometry for a stationary charged black hole, illustrating that 
in order to fall through the inner horizon, ingoing (positive energy) and outgoing (negative energy) streams must fall through 
separate ingoing and outgoing inner horizons into causally separate regions of spacetime. 

the case of two fluids, outgoing charged baryons could stream relativistically through ingoing neutral dark matter, 
and mass inflation duly occurred. 

I offer the following conceptual picture of the physical mechanism that causes mass inflation inside black holes. 
First, in order to drop from the superluminally infalling region of the black hole through an inner horizon into a 
subluminal region, it is necessary that ingoing and outgoing fluids drop through separate ingoing and outgoing inner 
horizons into causally separated pieces of spacetime, as seen most clearly in a Penrose diagram, such as the Penrose 
diagram of the Reissner-Nordstrom geometry shown in Figure [TJ If the ingoing and outgoing subluminal spacetimes 
were causally connected, then time would go in one direction for the ingoing fluid, and in the opposite direction 
for the outgoing fluid, which cannot be. This is most evident in stationary or self-similar black hole spacetimes, 
where there is a Killing vector associated with (conformal) time-translation symmetry, and that Killing vector must 
point in one direction for the ingoing fluid (positive energy geodesies), and in the opposite direction for the outgoing 
fluid (negative energy geodesies). The only way for the ingoing and outgoing fluids to pass into causally separated 
subluminal spacetimes is to move faster than the speed of light through each other, which is impossible. In their 
attempt to exceed the speed of light, the ingoing and outgoing fluids counter-stream ever faster, and in so doing they 
create a growing positive radial pressure. Mass inflation begins when the gravitational force produced by the radial 
pressure of the counter-streaming fluids exceeds the gravitational force produced by the mass of the black hole. The 
gravitational force produced by the radial pressure is inwards, but the inward direction is in opposite directions for 
the ingoing and outgoing fluids, being towards the black hole for the ingoing fluid, and away from the black hole 
for the outgoing fluid. The result is an exponential feedback instability, in which the increasing counter-streaming 
velocity increases the pressure, which increases the gravity, which accelerates the fluids ever faster through each other. 
Mass inflation thus appears to depend on an intriguing alternation between gravitational attraction and repulsion: 
first, there must be sufficient gravitational attraction to produce gravitational collapse and an event horizon; next, 
there must be sufficient gravitational repulsion to bring the black hole to the brink of having an inner horizon; and 
finally, mass inflation produces an exponentially growing inward gravitational attraction sourced by the pressure of 
counter-streaming ingoing and outgoing streams. 

Most published analytic and numerical investigations of mass inflation have considered feeding a spherical charged 
black hole with a massless scalar field, usually uncharged [1, Si, [11 [ll|, [11, III [13, [11, [^, [11, |4l|, ill, [s^ but sometimes 
charged jH, [2^, [H, J^, i47|, i7^, [87] . These investigations have confirmed that mass inflation occurs generically when a 
spherical charged black hole is fed with a massless scalar field. Since a massless scalar field supports waves moving at 



3 



the speed of light, it is consistent to hypothesize that mass inflation in the scalar field arises physically from relativistic 
counter-streaming between ingoing and outgoing waves. This idea provides one of the principal motivations behind 
the present paper: can one see explicitly from a perturbation analysis how counter-streaming between waves moving 
at the speed of light near an inner horizon might account for mass inflation? The present paper does not answer this 
question, but it lays a foundation for being able to address the question in future work. 

Gravitational waves are of particular relevance to the issue of mass inflation, because even if there were no other 
way in which relativistic counter-streaming could occur (perhaps there is no massless scalar field, and perhaps the 
fluid is so opaque that electromagnetic waves cannot propagate), certainly gravitational waves will allow it. However, 
spherically symmetric gravitational waves do not exist - they must be of at least quadrupole (l = 2) order - so it is 
necessary to go beyond spherical symmetry to admit gravitational waves. Of course, as remarked for example by [6^ . 
perturbations of a massless scalar field can be considered as a model of gravitational perturbations. However, this is 
something that, like the assumption that charge models rotation, should be checked. 

The purpose of this paper is to present the foundations of the perturbation theory of spherically symmetric, self- 
similar black holes. The work presented here is a natural generalization of the extensive literature on perturbations of 
spherically symmetric stationary black holes initiated by Regge & Wheeler (1957) [82*1; see especially the monograph 
by Chandrasekhar (1983) and the reviews by [13, El, [91| and references therein; and more recently [13, Hj, [S, 
[ssl . [6^ . [tqI . [ssj . General relativistic spherically symmetric self-similar solutions were first applied to the problem of 
gravitational collapse by Ori and Piran (1990) [iq, following earlier work on general relativistic self-similar solutions 
in cosmology. Since then there have been many investigations; see especially the reviews 15, 17, 43, 54] and references 
therein; and more recently [1, H HI [tiI HH, [93| . 

There are two reasons why it is useful to consider self-similar solutions and their perturbations as a way to understand 
mass inflation better. The first reason is that, although the geometry in the vicinity of an accreting black hole may 
be well approximated by an empty stationary solution (the Reissner-Nordstrom geometry, in the case of charged 
spherical black holes), the geometry inevitably departs hugely from the empty stationary geometry in the region near 
the inner horizon where mass inflation takes place. The simplest way to model this departure from the stationary 
case in a fully non-linear fashion is with self-similar solutions. The second reason is that numerical gravitational 
simulations are notoriously plagued with numerical instabilities (l^ . [30l . [sil . [s^ . The mass inflation instability is a 
real instability, not a numerical one, and it is essential that numerics be under tight control to explore mass inflation 
reliably. Self-similar solutions offer a numerically robust and fast approach to modeling mass inflation with some 
degree of flexibility. Although much of the literature on self-similar black holes is restricted to the case of a single, 
ideal fluid [1^ [l3, [H, H^l, the scope of self-similar solutions is actually much broader than this. It is possible to 
consider non- ideal fluids, for example with a finite conductivity [1^, and multiple, possibly interacting, fluids [ioj . 
Thus self-similar solutions provide a broad stage upon which to carry out investigations of mass inflation. 

Several papers have considered spherical perturbations to spherical self-similar spacetimes; see the review by [43} . 
and more recently (ssl . [6l| . By comparison, there have been few investigations of non-spherical perturbations of 
spherical self-similar spacetimes. The investigations that have been done have been based on the Gerlach and Sengupta 
(1979) [32*1 framework, which treats perturbations of arbitrary spherically symmetric spacetimes, not necessarily 
stationary or self-similar (see also the independent work by [49| on arbitrary spherical spacetimes). Gundlach and 
Marti'n-Garcia [1^ review the Gerlach-Sengupta formalism, and use it to analyze perturbations of the Choptuik [2ll[58| 
critical solution that describes the discretely self-similar collapse of a scalar field at the threshold of gravitational 
collapse (see [S^l for a review of critical gravitational collapse). Gundlach and Marti'n-Garcia 38, 57J extend the 
Gerlach-Sengupta formalism to develop equations that could be used to investigate the spherical collapse of a star 
numerically. Frolov [3l| used the Gerlach-Sengupta formalism to investigate perturbations to the Roberts [11] self- 
similar solution, an analytic solution describing the gravitational collapse of an uncharged scalar field. Nolan et al. 
[tH used the Gerlach-Sengupta formalism in their study of perturbations of self-similar Vaidya spacetimes admitting 
naked singularities. The approach of the present paper is more restricted than that of Gerlach-Sengupta, but has the 
advantage that it yiel ds a complete separation of variables of the wave equations. 

Nolan et al. [69l.l70l.F7H consider spherically symmetric self-similar collapse models that produce naked singularities. 
In these models, gravitational collapse leads to the formation of a naked singularity, which subsequently becomes 
covered by an event horizon, as illustrated by the Penrose diagram in Fig. 2 of [70|. The naked singularity is a source 
of unpredictability, and the Cauchy horizon, defined in the usual way to be the boundary of predictability, is thus 
a null cone that extends outwards from the formation event of the naked singularity. [69l |70| show that the flux 
of (arbitrary, non-spherical) scalar radiation impinging on this Cauchy horizon is finite, and |7l| show that, if the 
background self-similar solution is that of null dust (Vaidya geometry), then the flux of gravitational and matter 
perturbations remains similarly finite on the Cauchy horizon. Nolan et al. conclude that the finite flux points to the 
stability of the Cauchy horizon in this situation. This conclusion is consistent with the conceptual picture described 
above. Mass inflation, as envisaged by Poisson & Israel, must occur within a superluminally infalling region beneath 
an event horizon, for only then can the inward direction (towards smaller radius) point in opposite directions for 
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ingoing and outgoing streams. Thus the Cauchy horizon considered by [g^, [tOi l7l[. which Hes outside the event 
horizon, is not expected to be subject to the mass inflation instabihty. Curiously, [71| go on to find that the flux 
diverges on the event horizon that eventually covers the naked singularity. The physical origin of this divergence is 
not clear. 

The most fundamental result of this paper is a generalization of the Tcukolsky spin-s wave equation 0, [ssl . [soj 
to the case of spherically symmetric, self-similar black holes. There is no additional restriction on the form 
of the unperturbed energy-momentum tensor beyond spherical symmetry and self-similarity. As in stationary 
black holes, eigenfunctions of the wave equations for the propagating components of spin-s waves are separable, 
taking the form [(1 + V)/{1 — y)]^"^^ A^'^''+^^e^*'^*s'(/'w/(r) sYim{0, (jy) in conformal polar coordinates {t, r, 6, 0}, where 
[(1 + V)/{1 — 1/)]^*^^ is a radial Lorentz boost factor (which equals unity in a frame at rest in the conformal frame), 
A is a conformal factor, a; is a conformal frequency, are spin-s spherical harmonics, and s'>Puji{r) are radial 

eigenfunctions of the following generalization of the Teukolsky equation 

|- ^ + (^'^ ± + [Kl + 1) - - (1 + s^){R - P) - 2\s\Pi_ + 2(1 + 2s2)M] a| = sources (1) 

where the ± sign in front of k is -f for ingoing, — for outgoing waves. Here r* is a generalization of the Regge- 
Wheeler [s^l radial coordinate; A is the dimensionless horizon function; k = ^dln A/ dr* ; R = AnA^p, P = AnA^p, 
and P± = AttA^px^ are the dimensionless energy density, radial pressure, and transverse pressure; and M = m/A is 
the dimensionless interior mass, with m the interior, or Misner-Sharpe [59|, mass. The present paper considers only 
integral spin waves - gravitational waves (spin 2), electromagnetic waves (spin 1), and scalar waves (spin 0) - so it is 
not guaranteed that equation ([1]) holds for half-integral spin, although it may well do so. 

It should be emphasized that in the general self-similar case of non-vanishing background energy-momentum, the 
system of wave equations is inextricably coupled, so that waves of one kind excite perturbations which in turn provide 
a source for waves of another kind. Thus as long as waves are present, the source terms on the right hand side of 
equation ([1]) cannot in general be made to vanish, which is unlike the case of the Schwarzschild or Kerr geometries. 

The structure of this paper is as follows. Section [IT] describes the unperturbed self-similar background, ijllll sets 
up the Newman-Penrose formalism [H, [63] , and jjlVI uses the Newman- Penrose formalism to develop the equations 
governing perturbations of the self-similar background. Sections |Vl IVH and IVIII proceed to derive wave equations 
for gravitational, electromagnetic, and scalar waves. Section IVlIII summarizes the results. Because of the length and 
technicality of the paper, application to problems of interest, notably mass inflation near the inner horizon, is deferred 
to subsequent work. 

The conventions of this paper are for the most part those of Misner, Thorne & Wheeler [63| (with the exception 
that the indices on the electromagnetic field tensor Fmn are swapped compared to MTW). All complicated equations 
have been calculated and checked with the algebraic manipulation program Mathematica. It is to be hoped that if 
there are errors in the equations, then they are errors of transcription, and they do not propagate. 

In many places in this paper, starting from ^III C[ symmetrically related pairs of equations are given as single 
equations with two sets of alternative indices. The notation is used partly for compactness, and partly to manifest 
the u <-^- w, -|- <-> — symmetry of the equations in the Newman-Penrose formalism 33]. To reduce the possibility of 
error, the algebra for each of a pair of symmetrically related equations was calculated separately, the pair were then 
merged, and finally the indices were checked for the expected symmetry. 

II. UNPERTURBED BACKGROUND 

This section presents the equations governing the spherically symmetric self-similar black hole solutions that provide 
the unperturbed background to the subsequent perturbation theory. In basic respects the development is similar to 
that of [3^, 113] , but the gauge choices that permit separation of variables in the wave equations prove to be different 
from the gauge choices of fssl, S3i so it is necessary to re-develop the formalism ab initio. The most fundamental 
difference from [s^ 14^] is to make the conformal character of the self-similar solutions explicit from the outset. 

A. Self-similarity and the choice of coordinates and tetrad 

The assumption of self-similarity is the assumption that the system possesses a conformal time-translation symmetry. 
This implies that there exists a conformal time coordinate t such that the geometry at any one time is conformally 
related to the geometry at any other time 



ds^ = a{tfgl^}dx>'dx' 



(2) 
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in which the conformal factor a{t) depends only the conformal time t, and the conformal metric coefficients g^^} are 
all independent of the conformal time t. The condition ^ on the metric is a necessary but not sufficient condition 
for self-similarity: self-similarity requires that the entire system, including notably the energy-momentum tensor, 
possesses conformal time-translation symmetry. For example, the Friedmann-Robertson- Walker metric of cosmology 
satisfies equation ([2]), but in general does not evolve in a self-similar fashion. 

The assumption of spherical symmetry implies that the metric can be reduced further to the form, in conformal 
polar coordinates x'^ = {t,r,9,(j)}, 

ds^ ^a{tf \g[f {r)dt^ + 2 g[$ {r)dtdr + gi") {r)dr^ + r^do^l (3) 
in which do^ = dO^ + sin^9 dcjP' is the metric of the surface of a unit sphere, and the conformal metric coefficients g[f , 

(c) (c) . 

glr , and gW are functions only of the conformal radial coordinate r, independent of either the conformal time t or 
the angular coordinates, the polar angle 9 and the azimuthal angle (p. 

In place of the conformal factor a it is convenient to work with the conformal factor A 

A EE a{t)r (4) 

since it proves to be A^ not a, that appears ubiquitously in formulae. The fact that the angular part of the metric ([3]) 
is a^r^do^ = A^do^ implies that the conformal factor A is equal to the proper circumferential radius, with the defining 
property that the proper circumference measured by observers at fixed conformal time t is 2ttA. Being a proper 
quantity, A is coordinate (and tetrad) gauge- invariant, whereas neither the conformal time t nor the conformal radius 
r are themselves coordinate gauge-invariant. The conformal factor ^4 is a dimensional quantity, with units of length. 
In self-similar solutions, all quantities are proportional to some power of A, and that power can be determined by 
dimensional analysis. 

In this paper I choose to work in the tetrad formalism. Orthonormal tetrads are considered in this section, i jlll and 
associated Newman-Penrose [H, 113] null tetrads will be considered starting from mil\ 

Let denote the basis of tangent vectors of the coordinates. By definition, the scalar products of the tangent 
vectors constitute the metric 

9^^ gu= 5^1/ ■ (5) 

In the tetrad formalism, a set of frames is erected at each point of the spacetime, with axes 7„i7 the tetrad, whose 
dot products form the tetrad metric 7mn 

(6) 

In an orthonormal tetrad, as considered throughout this section, the axes are chosen to be locally inertial, so that the 
tetrad metric 7„i„ is the Minkowski metric rimn 

In a Newman-Penrose tetrad, considered in jjllll and thereafter, the tetrad metric is not Minkowski, but the tetrad 
metric components 7mn, equation (|53p . are still constants, so that for example equation (|30p below, and thus also 
equation (|33p . remains valid. In this paper dummy latin indices signify tetrad frames, while dummy greek indices 
signify coordinate frames. The axes 7„i of the locally inertial frames are related to the tangent vectors by the 
vierbein and its inverse e™^ 

— Q ^ ; — ^ ■ (8) 

The vierbein provide the means of transforming between the tetrad components and coordinate components of any 
4-vector or tensor object. For example, the tetrad components and the coordinate components of the 4-vcctor 
p = p^gf_i — p™jm are related by p™ — e^^p^. Tetrad components are raised, lowered, and contracted with the 
tetrad metric jmm while coordinate components are raised, lowered, and contracted with the coordinate metric g^jy. 

The most general vierbein that is manifestly spherically symmetric contains 6 free parameters [s^ . of which 2 
correspond to physically distinct spacetimes, 2 are associated with coordinate gauge freedoms of time and radius, and 
2 are associate with the tetrad gauge freedoms of a Lorentz boost in the radial direction, and a spatial rotation about 
the radial direction. I choose to eliminate one of the coordinate gauge freedoms by choosing the conformal factor A, as 
above, equation ([?]), to be the proper circumferential radius, and I choose to eliminate the tetrad freedom associated 
with spatial rotations about the radial direction by aligning the tetrad axes with a non-spinning locally inertial frame 
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in radial free-fall. Specifically, I choose the vierbein to take the following form, whose slightly complicated form yields 
a somewhat simpler inverse vierbein, equation (jlOp . 



1 
A 



( 



V 

















— 





1/sin^ 



(9) 



The choice of vierbein in equation ^ is more general than necessary, but it is convenient to defer the choice of gauge 
of the conformal radius r (which will in mil Al be dictated by the requirement that eigensolutions of wave equations 
be separable in the coordinate r), and to retain the freedom of an arbitrary Lorentz boost in the radial direction (so 
that, for example, one could choose the tetrad to be comoving with the center-of-mass frame of infalling matter). The 
vierbein ^ contain 4 free parameters, which in self-similar solutions turn out to constitute two dimensionless tetrad- 
frame 4-vectors ^™ = {C*,^'^jO,0} and A'" = {A*, A'', 0, 0}. Dimensionless means a function only of the conformal 
radius r, independent of the conformal factor A. The conformal factor A itself is not a free parameter, being fixed 
equal to the circumferential radius, equation ([4]). The quantity ^ • A = — ^*A* -I- ^''A'' in equation ([9]) is the scalar 
product of the dimensionless 4-vectors ^ and A. The inverse of the vierbein ^ is 



= A 



/C* A7r 

C A7r 

10 

V sin( 



(10) 



The metric ds^ = jmnG™ f^e"^ i.dx'^ dx"^ corresponding to the inverse vierbein e™^ of equation PU)) is 



ds' - A' 



^^dt+—drY 
r 



-dr + Cdt 



do'' 



(11) 



and it apparent that the assumption of self-similarity, equation is equivalent to the hypothesis that, as claimed 
above, ^™ and A™ are dimensionless, independent of the conformal time t. The determinant of the metric is equal to 
the square (1/e)^ of the determinant 1/e of the inverse vierbein, which is 



i=A4^sin( 
e r 



Directed derivatives dm are defined to be the spacetime derivatives along the axes 7m of the tetrad frame: 



d„ 



7m ■ d 



(12) 



(13) 



where d = g'^d/dx'^ = g'^'^g,yd/dx'^ is the invariant spacetime vector derivative. The directed derivatives dm depend 
only on the choice of tetrad frame, and are independent of the choice of coordinate system. Unlike the coordinate 
partial derivatives d/dx^, the directed derivatives dm do not commute, equation (132]) below. 
The coordinate 4- velocity of an object at rest in the tetrad frame is 



= dtx" = e7 = 



A^- X 



(14) 



If the tetrad frame were chosen to be at rest in the similarity frame, so that dtr = 0, then this would correspond to 
the gauge choice 7 = 0- However, I choose to retain the freedom of an arbitrary radial Lorentz boost, that is, I allow 
7 to be arbitrary. 



B. Gauge-invariant scalars 



The fact that quantities in self-similar solutions are, modulo powers of the conformal factor A, independent of 
conformal time t (with vanishing conformal time derivative d/dt) means that the tangent vector gt is a conformal 
Killing vector: 

..•a = |. (15) 
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The tetrad- frame components of the Kilhng vector are A^™, because 

7" • 9t = e"* = AC (16) 

or equivalently 

^ - Ardra ■ (17) 

This explains the choice of the quantities i^™ in the definition ([9]) of the vicrbein, and also demonstrates that is a 
tetrad-frame 4-vector, as claimed. 

Minus the square of the dimensionless Killing vector ^™ defines the horizon function A, a coordinate and tetrad 
gauge-invariant dimensionless scalar: 

A = -e^r = -A-^gu . (18) 

According to the metric ((TT|l . the proper time r experienced by observers at fixed conformal position in the black hole, 
dr — d6 = d(f> = 0, is related to the conformal time i by dr = AA^^^dt, the vanishing of which defines the positions of 
horizons. The horizon function A is positive outside the horizon, zero at the horizon, and negative inside the horizon, 
reflecting the fact that the Killing vector gt is timelike outside the horizon, null at the horizon, and spacelike inside 
the horizon. 

Besides the dimensionless Killing vector , there is another coordinate gauge-invariant dimensionless tetrad- frame 
4-vector that arises naturally in self-similar solutions, the radial 4-gradient (3rn = {Pt, f3r, 0,0} defined by 

Pm = dmA . (19) 

Physically, Pt and Pr are the proper time and radial derivatives of the circumferential radius A as measured by 
observers at rest in the tetrad frame. That Pm is coordinate gauge-invariant follows from the fact that it is a tetrad- 
frame derivative of the coordinate gauge-invariant quantity A. The fact that P„i is dimensionless follows from the 
fact that, if f{r) is some arbitrary dimensionless function, then 

d,nAf{r) = Ad„J{r) + Pmf{r) (20) 

is a sum of two terms, the first of which is manifestly dimensionless because Admf{r) = Acm^df /dr and Cm'^ oc l/A, 
equation ([9]), from which it follows that the second term Pmf{i") must also be dimensionless, so that Pm itself must 
be dimensionless as asserted. Expressions involving directed derivatives of powers of the conformal factor A times 
dimensionless quantities, such as on the left hand side of equation (j20p . appear ubiquitously in the theory of self-similar 
black holes, for example in the definitions of the tetrad-frame connections and Riemann tensor. 

As an aside, it is worth mentioning that the argument from equation (|20p that Pm is dimensionless fails for the 
Friedmann-Robertson- Walker metric, which possesses spatial translation symmetry, so that physical quantities depend 
only on conformal time t, not on position, and consequently /(r) in equation (j20p cannot be taken to be arbitrary. 
Thus in the Friedmann-Robertson- Walker metric a more general class of self-similar solutions exists. However, the 
present paper concerns itself with self-similar solutions of black holes, not of cosmology, so for this paper Pm is indeed 
dimensionless. 

The scalar product £,"^Pm of the Killing vector ^™ with the radial 4-gradient defines a second coordinate and 
tetrad gauge-invariant dimensionless scalar v 

Since a(t) is a function only of conformal time i, equation ([2]), the last expression of equation (j2ip . v = (daj dt) j a, 
shows that v is a function only of conformal time t. But the fact that v is dimensionless (a function only of r) in 
self-similar solutions implies that v must also be independent of conformal time t. It then follows that v must be a 
constant, independent of either t or r. In [s^ |43| v was set to unity, which corresponds to a certain gauge choice of 
the scaling of conformal time t. Here instead v is allowed to be an arbitrary constant, which has the virtue of making 
it transparent how to take the limit of a stationary black hole, which corresponds to v — > 0. Physically, the constant 
V is proportional to the expansion velocity of the black hole at some fixed (but arbitrary) conformal radius r, and 
thus a measure of the rate at which the black hole is accreting and growing. Specifically, the velocity v is related to 
the rate dA/ dr at which an observer at rest in the conformal frame measures the circumferential radius A to vary 
with proper time dr — AA^^^dt by 

y = A-i/2^. (22) 
dr 
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In terms of the expansion velocity v and the 4- vectors ^"^ and A™ in the vierbein ([9]), the components /3m = 
em^dA/dx^^ of the radial 4-gradient are 

A third coordinate and tetrad gauge-invariant dimensionless scalar, besides the horizon function A and the expansion 
velocity v, is given by the square /3™/3m of the radial 4-gradient In terms of this a gauge-invariant dimensionless 
scalar, the dimensionless interior mass M can be defined by 

M = i (1 - PraD ■ (24) 

The dimensionless interior mass M is related to the the interior, or Misner-Sharpe [H^, mass m by M = ra/A (this 
notation differs from that of [s^ l40l|. who defined M to be the interior mass, as opposed to a dimensionless version 
thereof). The physical interpretation of to = AM as an interior mass emerges from the spherically symmetric Einstein 
equations. 



C. Tetrad-frame connections 

In the tetrad formalism, the covariant derivatives -D„ of the covariant, pk = ■ P, and contravariant, p'' = ■y'^ ■ p, 
tetrad-frame components of a 4-vector p = p^g^ = p"^'fm are 

DjiPk = dnPk - T^TnPm 

Dnp" = dnp'^+Tl^p"^ (25) 

where the tetrad- frame connection coefficients r^„, also known as the Ricci rotation coefficients, are defined by the 
directed derivatives (9„ of the tetrad axes 'jm 

dnlra = T^nlk (26) 

in much the same way that the usual coordinate-frame connection coefficients, the Christoffel symbols TJ^^^, are defined 
by the coordinate derivatives d/di, of the coordinate tangent vectors 

^ - . (27) 

Let denote the vierbein derivatives 

(7fe ^e'^de'^-e'^e'' ^^"'^ (28) 

The fact that 7^ = Cm^g^ implies that the tetrad-frame connection coefficients defined by equation (pS)) are 
related to the coordinate-frame connection coefficients F^^ defined by equation (P7)) by 

mn ~ ^mn ~ K^m c„ -L . \^-') 

The definition (|26p and the fact that the tetrad metric coefficients are constant, 9„(7fc • jm) = dnjkm = 0, implies 
that the tetrad-frame connection coefficients with all indices lowered Tkmn = Iki^mn are antisymmetric in their first 
two indices. 



kmn 



— —^mkn (30) 



which expresses mathematically the fact that Tkmn for each given final index n is the generator of a Lorentz 
transformation between tetrad frames parallel transported along the axis 7„. The tangent vectors can be regarded 
as coordinate derivatives of the invariant 4-vector interval dx = g^j^dx^, that is, g^ — dx/dx^, and the commutativity 
of partial derivatives, dg^/dx'^ = d^x/dx'^dx^ = d^xjdx^dx^ = dg^/dx^, implies the usual no-torsion condition of 
general relativity that the coordinate-frame connections are symmetric in their last two indices 

r^i/ = Tiz/j • (31) 
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Unlike coordinate derivatives, the directed derivatives do not commute. The commutator of directed derivatives is 

[dm dm] = (r,„„ — r„„j)9fc = ((i„„j — <^rim)^fc (32) 



the middle expression of which follows from equation (|26p . and the last expression from equation (|13p . From the 
relation the Lorentz antisymmetry (PU]) . and the no-torsion condition ([?!]) . it can be deduced that the tetrad- 
frame connection coefficients Tkmn can be expressed in terms of the vicrbcin derivatives dkmn = Ikidlmm equation (j28p . 
as 



fcmn 



which provides an explicit way to compute Tkmn from the vierbein Cm^- 

For the vierbein e™'^ given by equation the non- vanishing tetrad- frame connection coefficients Tkmn are 

= /it/A , (34a) 

Ft.t = , (34b) 

r^t^ = Tete = /?tM , (34c) 

T^r4, = Te.e = Pr/A , (34d) 

F^e^ = cote/ A , (34e) 



where 



' C • A V 5* dr J ' ^-Xy dt dr 

Physically, the connection coefficient ht/A = Ttrr is minus the local "Hubble parameter" of the radial flow, that is, 
minus the proper radial gradient, —D^vT, of the proper radial velocity between objects each of which is comoving 
with the tetrad frame at its position; while hr/A = Ttrt is the proper gravitational force, or equivalently minus the 
proper acceleration, —Dtu^, experienced by an observer at rest in the tetrad frame. These assertions follow from the 
fact that observers at rest in the tetrad frame have 4-velocity = {1, 0, 0, 0}, so that 

Dru'' = dru'' + T^y = F[^ = ~ht/A , (36a) 
Dtu"^ = dtu'' + FJ-^u* = Tl = -hr/A . (36b) 



D. Riemann, Ricci, Einstein, and Weyl tensors 



The Riemann tensor Rkimn in the tetrad frame is defined in the usual way by the commutator of the covariant 
derivative, RkimnP" = [Dk, Di]Pm, and is given in terms of the tetrad- frame connection coefficients by 



R 



klmn 



d,T 



nmk 



dkT 



k^ nml 



yj r 

m/-^ jnk 



y] r 

mk^ jnl 



(F^- 



Ik 



ri.)r„. 



(37) 



which has two extra terms (the last two) compared to the usual coordinate expression for the Riemann tensor in 
terms of Christoffel symbols. The Ricci tensor and scalar are given by the usual contractions of the Riemann tensor, 
Rkm = l''^^Rkimn and R = ^''"^Rkm, and the Einstein tensor is given by the usual expression in terms of the Ricci 
tensor and scalar, Gkm — Rkm ~ \Rlkm- The Weyl tensor Ckimn is defined in the usual way to be the totally 
trace-free, or tidal, part of the Riemann tensor Rumm 



Cklmn = Rklmn — \ [ikmRln — IknRlm + JlnRkm — llmRkn) + \ {jkmlln " Iknllnx) 

In the present case, define the symbols F, R (not the Ricci scalar!), P, and Pi by 

F = Adtfir + PtK = Adrf3t + I3rht , 

R = -Adrf3r- Ptht + M = -^[dr{AM)- (3tF] , 



p 

p± 



-AdtPt - Prhr 



M 



1 



[dt{AM) + PrF] 



1 



A{drP + dtF) -hr{R + P)- 2htF 



1 r 

Wt 



A{dtR + drF) -ht{R + P)- 2hrF 



(38) 

(39a) 
(39b) 

(39c) 

(39d) 



10 



The definitions ((39)) of F, R, P, and P± differ from the corresponding definitions in [s^illol in that appropriate factors 
of the conformal factor A have been included so as to make F, R, P, and P± dimensionless, that is, functions only of 
the conformal radial coordinate r in self-similar solutions. In terms of these quantities, the non- vanishing components 
of the Einstein tensor Gmn are 

Gtt = 2R/A^ , (40a) 
Gtr - -2F/A^ , (40b) 
Grr = 2P/A^ , (40c) 
G^^ = Gee - 2PJA^ . (40d) 

The Einstein equations 

Gmn = SnTmn (41) 

imply that the dimensionless quantities R, F, P, and P± are related to the tetrad-frame energy density p = T", 
radial energy flux / = r**", radial pressure p = T^^ , and transverse pressure p± = T^^ — T"^"^, by 

R = 'iirA^p , F = AttA^J , P = AnA^p , P±_ = AtiA^p±_ . (42) 

The non-vanishing components of the Weyl tensor Gkimn are 

Gtrtr = —Ce,pe<i) — —2Gtete = —2Ct,pt<p — 2Grere — 2Gr4,r4, — C/A^ (43) 

where C is the dimensionless Weyl scalar 

C=l{R-P + Pj_-3M) . (44) 



III. NEWMAN-PENROSE FORMALISM 



The perturbation theory presented in i jl VI follows the Newman-Penrose formalism 0, [U [T^I , and this section sets 
up the necessary apparatus. Subsection llll Al dcfines the background radial gauge with respect to which wave equations 
prove to be separable. Subsection llll Bl introduccs the Newman-Penrose tetrad. Subsections IIII C l lIII El define angular 
and radial raising and lowering operators, in terms of which wave operators turn out to be most elegantly expressed. 
Finally, ijlllFI discusses the asymptotic behavior of wave solutions in a hypothetical asymptotically flat spacetime at 
large distance from the black hole, a question that can be addressed independently of the black hole itself. 



A. Radial gauge and the Regge- Wheeler coordinate 



Up to this point the gauge of the conformal radial coordinate r in the unperturbed background has been left 
unspecified. It turns out that the requirement that eigenfunctions of wave equations be separable functions of the 
conformal coordinates t, r, 0, and (j) forces the radial coordinate r into a certain gauge. It should be commented that 
in much of the rest of this paper the background radial gauge is not actually fixed until the last moment. That is, the 
definitions (|56aP and (|69p of the radial operators Vv and s'Dv are independent of the background radial gauge: only 
when these operators are converted to coordinates, equation (|57ap . is the radial gauge actually fixed. 

Recall that the coordinate time derivative is the directed derivative along the Killing vector A£^"^ = A {^*, 0, 0} 

^ = ^r5™. (45) 

Let X™ = {S.^, £,* ,0, 0} denote the 4- vector in the r-t plane orthogonal to the Killing vector, with the same magnitude 
as The directed derivative along this orthogonal direction is 

Separability of the wave equations turns out to require that the directed derivative along the orthogonal direction 
X™, equation (^5]) . be a purely radial coordinate derivative. This imposes the gauge condition that ^*Ar + ^''A^ = 0. 
The corresponding values of A™ are 

A™ = ^ — . (47) 

[(1 -2M)A + v2]i/2 
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In this gauge, the directed derivative in the direction orthogonal to the KiUing vector ^"^ is a purely radial 
derivative, as desired. 



The coordinate metric in this gauge reduces to the diagonal form 

A 



21 1/2 



d_ 

dr 



Adt^ 



[(1 - 2M)A + v2] r-2 



dr"^ + do^ 



(48) 



(49) 



Note that even though the coordinate metric g^i, is diagonal in this gauge, the vierbein Cm'^ is not diagonal, thanks 
to the (desirable) freedom of an arbitrary Lorentz boost in the radial direction. 

In self-similar evolution, the quantities A, v, and M, equations ([T5)) . (ED), and ([M]) . are all dimensionless, which 
is to say functions only of the conformal radius r, and a Regge- Wheeler ^] (tortoise) coordinate r* can be defined 
such that 



_d_ 

dr* 



In terms of the Regge- Wheeler coordinate r* , the metric then simplifies further to 

ds^ = A^A{-dt^ + dr*^)+do^] . 



(50) 



(51) 



Stationary solutions are obtained in the limit of zero expansion velocity, v — s- 0. In this case the time-dependent 
part a{t) of the conformal factor A = a{t)r, equation ([4]), can be set to unity without loss of generality, a(t) = 1, so 
that it is consistent to take A — r in stationary solutions. 

The Reissner-Nordstrom geometry is not only stationary but also empty aside from a central singularity and a 
static electric field. To convert to the Reissner-Nordstrom geometry, set A = (1 — 2M)/r'^, the dimensionless interior 
mass M being related to the black hole's constant mass m, and charge q, by M = m,/r — jq^/r^ ■ The conformal 
time t and conformal radius r then coincide with the usual Reissner-Nordstrom time t and circumferential radius r. 



B. Newman-Penrose tetrad 



Newman-Penrose null tetrads [33, [63] are particularly well adapted to exploring fields that propagate at the speed 
of light [77| . Define the Newman-Penrose null tetrad {7^, 7„, 7+, 7_} in terms of the orthonormal tetrad {7^, ^r,le,lct>\ 

by 



The Newman-Penrose tetrad metric is 



Iv 
lu 

7+ 
7- 



^(7.+>), 
^(7.->), 

;^(7e -i70) ■ 



7m • In 



/ -10 
-10 
1 
V 10 



(52a) 
(52b) 
(52c) 
(52d) 



(53) 



with indices to, n running over d, u, -)-, — . 

It is useful to define the operation of angular conjugation, designated * (to be distinguished from complex 
conjugation *), as flipping the azimuthal tetrad axis: 



70 = -!<!> 



(54) 
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In a Newman-Penrose tetrad, angular conjugation corresponds to swapping the angular axes, +<-!■—, 

ll = lv , lu = lu , ll=l- , 7- = 7+ • (55) 

If a quantity is real in an orthonormal tetrad, then in a Newman-Penrose tetrad its angular conjugate is the same as 
its complex conjugate, but with waves it is convenient to allow a complex time dependence ~ g-*'^* (and an associated 
complex radial dependence), in which case angular conjugation is not the same as complex conjugation. The effect 
of angular conjugation is to swap covariant angular indices + <-> — on tensors (such as the Riemann tensor) and 
other natively real objects in a Newman-Penrose tetrad (such as the tetrad connections). In general, an object can 
be resolved into polar (p) and axial (a) parts defined by the property that the polar and axial parts respectively do 
not and do change sign under angular conjugation. For objects that are real in an orthonormal tetrad, the polar and 
axial parts correspond to the real and imaginary parts of the object in a Newman-Penrose tetrad. 



C. Commuting radial and angular operators 

The radial derivative operators Vy and I?„ and angular operators X>+ and I?_ defined as follows, equations (|56p. 
play a fundamental role because they provide the basis for the definitions of radial and angular raising and lowering 
operators [equations ([69)1 and ((6T|l below] which appear ubiquitously in the equations describing perturbations of 
spherically symmetric self-similar black holes: 

Vv = \dvA^ , (56a) 
X>± = Ad^ . (56b) 
In the radial gauge specified in ijlll Al equation ((50|) . the operators I?™ become 

The radial operators "D^ and I?„ commute with the angular operators I?+ and 'D_ (whereas the radial directed 
derivatives and i9„ do not commute with the angular directed derivatives 9+ and d_ ) . The placement of factors of 
A in equation (|56ap for Vv is such as to eliminate factors of A in the relation ([70)1 below. 

The eigenfunctions of wave equations are mutually orthogonal. The measure with respect to which inner products 
of eigenfunctions are defined is the invariant 4-volume element dl^x ~ e^^dtd(f>d9dr, where e is the determinant of the 
vierbein Cm^. In the present case the invariant 4-volume d'^x is a product of radial and angular parts do and do 

d'^x = do do (58) 

given by 

do = A^Adtdr* , (59a) 

do EE sin 6 dOdcj) . (59b) 

The contravariant components of the operators given by equations ([5^ are I?" = — I?" — —V^, and I?+ = 
2?" = 2?+, in accordance with the Newman-Penrose metric (j53p . and the Hermitian conjugates of these with respect 
to the radial and angular measures given by equations (|59p are 

pti = d^^A = —V^smO. (60b) 

sin 6 



D. Spin raising and lowering angular operators, and spin-weighted spherical harmonics 

This subsection summarizes the familiar properties [35| of the spin raising and lower angular operators and the 
spin- weighted spherical harmonics introduced by Newman and Penrose [67j . The aim is not to introduce any new 
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ideas, but rather to establish notation, and because the pattern provides a template for the spin radial operators 
presented in ^III El below. 

Define spin raising and lowering ang ular operators sl?+ and (these are, modulo a factor of — 1/V2, the same 
as Newman and Penrose's S and § [35l Isst) by 

sV+ = (sin 6I)''X>+ (sine)-" , (61a) 
sX>_ = (sin6l)"''X>_ (sin6l)" . (61b) 

The Hermitian conjugates of the contravariant components of these operators are 

^ _ ^ (62a) 

sV^- - -s-iV+ . (62b) 

The products of the raising and lowering operators s'D± with their contravariant Hermitian conjugates ^2?^* are the 
Hermitian operators 

2,2?t+,2?+ = ,L''-s{s+\) (63a) 
2sV^- sV_ = sL^-s{s-l) (63b) 

where sL'^ is the squared total angular momentum operator 

1 d 1 

= _ sin 61— + — ^ {m^+s^+2mscos9) (64) 

sm oO oO sm 6 

with d/d(j) = iJTi. 

An important feature of the operators s^?^^ s2^+ and s^?^" s^?- is that they differ by a constant, equations ((63|) . 

^2?^ „ ^pt- = _5 ^ (65) 

since it is this property that ensures that the ladder of raising and lowering angular operators s'C'± yields a mutually 
consistent set of eigenf unctions. The operator ^1?^+ represents the operation of raising followed by lowering, 
and thus connects spin-s with spin-(s+l) functions, while ^I?^" s'D^ represents the operation of lowering followed by 
raising, and thus connects functions on the adjacent level, spin-s with spin-(s— 1). The fact that the two operators 
differ by a constant, equation ([65| . implies that the spin-s eigenfunctions are the same for both levels. This is true for 
any spin s, and so it follows that, as asserted, the ladder of raising and lowering angular operators yields a mutually 
consistent set of eigenfunctions. 

The eigenfunctions of the squared total angular momentum operator sL"^ are the spin- weighted spherical harmonics 
sYiyn{d, (/)), with I > |s|, |m|, satisfying 

,Yira{0, <P) = l{l + 1) sYun{e, 0) . (66) 

The spin s, harmonic number I, and azimuthal number m must be either all integral, or all half-integral, although 
only integral spin is considered in this paper. The spin raising and lowering angular operators raise and lower by one 
unit the spin s of the spin-weighted spherical harmonics 



{lTs){l±s^\)Y'^ 

s±lYira{0,4>) . (67) 



The normalization factor in equation ([S7)) follows, modulo a conventional choice of sign (per [la, [6l|), from the fact 
that the raising and lowering operators sl?+ and —g+iD^ are Hermitian conjugates of each other, equation (|62p . and 
that their product —s+iD_ sl?+ has eigenvalue (Z— s)(Z+s+l)/2, in accordance with equation (|63ap . Spin-harmonics 
of opposite spins s are related by 

sYi„, = (-)'"+^_.y,*_^ . (68) 



E. Spin raising and lowering radial operators 

As will be seen in mVl radial operators and s^'u analogous to the angular operators 3!?+ and of the 
previous subsection, mil Dl appear ubiquitously in the perturbation theory of spherically self-similar black holes. The 



14 



structure of the radial operators resembles that of the angular operators, and it is convenient to refer to the radial 
operators as spin raising and lowering operators. However, whereas the angular operators yield a consistent ladder 
of angular eigenfunctions, the radial operators do not, because whereas in the angular case the Hermitian operators 
connecting adjacent levels of spin differ by a constant, equation (|65p . in the radial case they do not, equation (j75p 
(except in the case that the Weyl scalar C is constant, which occurs in flat space, where C = 0, as considered in the 
next subsection, §111 Fp . Thus the designation of s'Dv and s'Du as spin raising and lowering radial operators has a 
limited meaning. It is true that for each spin s, the Hermitian operator ^T)'^'" gDy has two eigenfunctions, and that s2?„ 
"raises" the "spin-s" eigenfunction to the "spin-(s+l)" eigenfunction, while sl^^" = s+i^u "lowers" the "spin-(s+l)" 
eigenfunction to the "spin-s" eigenfunction. But the eigenfunctions of operators g'^'^^ sT^v with different spins s are 
different, so there is no consistent ladder of spin radial eigenfunctions (except in the case that C is constant, considered 
in the next subsection, Fp . 

Define the spin raising and lowering radial operators gT^v and g'Du by 



gV, = icrv^ic)-' , 

The Hermitian conjugates of the contravariant components of these operators are 

gV^^ = g_,Vy . 



(69a) 
(69b) 

(70a) 
(70b) 



The products of the raising and lowering operators with their contravariant Hermitian conjugates s^^^ a-re the 
Hermitian operators 



2 s^v — 

2 qT)^ iiT^ii. 



1 



1 



d_ 

di 
d_ 

dt 



d_ 

dr* 
_d_ 

dr* 



c)lne" 

dr* 

(91ne 

dr* 



d_ 

dt 
d_ 
dt 



d_ 

dr* 
_d_ 

dr* 



d\nC 

dr* 

dlnC 
dr* 



These Hermitian operators can be re-expressed as 



2gV^W, = 



1 



1 



-s/2 



-s/2 



d^ 
d^ 



where d/dt = — icj, the dimensionless scalar k is 

1 9A 

'^^2Kd^ 

and C is the dimensionless unperturbed Weyl scalar, equation ((4^ 
in equations (|72p comes from 



+ [iuj + snf + s(l + QC)A 
+ {iuj + skY - s{l + QC)A 



A' 



A' 



%-\ A' 

i_ 



s/2 

^ 1 



d\nA 

dr* 



(71a) 
(71b) 

(72a) 
(72b) 

(73) 



The factor l-f 6C in the potential s(l -I- 6C)A 



1 dn 
A dr* 



= 1 + 6C 



(74) 



The Hermitian operators gV'^'" gVy and ^I?^" g'Du differ by a function, according to equations ([7^ . 

^2?^ _ ^pt« ^25^ = = s{l + QC) . 



(75) 

This difference, equation fTS]) . is generally not constant because the Weyl scalar C is generally not constant, and it 
follows that, in contrast to the angular raising and lowering operators g'D±^ the ladder of radial raising and lowering 
operators gX>v does not yield a mutually consistent set of radial eigenfunctions (except in the case that C is constant, 
such as the case C — Q considered in the next subsection, §111 Fp . 

In practice, the radial operators g'D'^'" gDy and s2?^" s'Du occur in the gravitational, electromagnetic, and scalar wave 
equations (fT72)) . (fT99)) . ([StTjl . ((269)) . and ((592)1 in combination with angular parts ^1?^+ and g'D^- g'D_ as follows, 
with positive spin s: 



1 



gV^^^gVu + ^gV^^^gV^^^^^^ 



dr* 



+ {iuj ± snf + + 1) - _ 6sC] A \ [ ^ 



±s/2 



A^ (76) 
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in which the total angular momentum operator s_L^ in the angular operators s'l^^'^ s1^+ and s1^^~ s2?-, equations (j63p . 
has been replaced by its eigenvalue 1(1 + 1), equation l|66p . The wave operators in the wave equations for gravitational, 
electromagnetic, and scalar waves prove to be the same as that given by equation (|76p . apart from the addition of 
further dimensionless terms, functions only of conformal radius r, in the factor inside square brackets multiplying the 
dimensionless horizon function A within the curly braces on the right hand side. 
The operator given by equation ([76|) has homogeneous solutions ^'±s 

(±,I?t- ^^Vn + is-D^^ is-D^) - (77) 

of definite conformal frequency uj and definite spherical harmonic numbers Im, separating as 

*±s = i^yC)^^/^A-' e-'"* ^Mr) isYimie, 0) (78) 

where ±s'0aj/(r) depends only on conformal radius r, not on conformal time t. The factors of (^^/^")Ts/2 ^j^g 
solutions \E'±s embody their dependence on radial Lorentz boosts. In a tetrad frame at rest in the similarity frame, 
where = and hence = the factor is just unity. In a tetrad frame moving at radial velocity V with respect 

to the similarity frame, so that i^' /^' = V, the factor (^"/^")Ts/2 jg simply the =Fs'th power of the special relativistic 
Doppler shift factor 



/^t) \ 1/2 


[(1 + ^)1 


1/2 


ley 







(79) 



Reduced spin raising and lowering operators sdv can be defined which operate only on the radial part si^uiiif) of 
the homogeneous solution 5*8, in view of the relation ((78)) . 

sdv = e*'^*A2(f7,f")('*±l)/^I?„ (^«/^«)-s/2^-2g-ic^t ^ (gQ) 

which simplifies to 

dv = i^) , sdv = A^'/^ dv AT"/2 . (81) 

" \/2A ^ dr* / " « 



F. Asymptotic wave solutions at large radius 

Self-similar solutions to black holes do not necessarily continue self-consistently to asymptotically flat, empty space 
at large radius [H, |4^. Nevertheless, it is convenient to suppose that the self-similar geometry does asymptote at 
large radius to a region where the Weyl scalar vanishes, C — 0, since this makes it straightforward to identify which 
solutions are to be identified as ingoing and which as outgoing. 

The solutions of equations ((73|) and ([74|) with C = 0, and with k and A tending to zero at large radius, are 

- = (82) 

In non-stationary self-similar solutions, where the expansion velocity is non-zero, v ^ 0, the Regge- Wheeler coordinate 
r* defined by equations (|48|) and (j50|) is related to the conformal radius r at large radius by vr* ^ Inr. This contrasts 
with the case of stationary solutions, where r* ^ r at large radius. 

Denote the radial part si^uii{f) of the solutions 5*5, equation (|78p . in the asymptotic region at large radius by sZi 

sipuji sZl at large radius . (83) 

The asymptotic solutions sZi are the homogeneous solutions of the operator inside braces in equation (|76p for the case 
where C — and k and A are given by equations (j82p : 



„ 2siuj 1(1 + 1) 



sZi^Q. (84) 



It evident from the form of equation (|84p that the solutions can be taken to be functions sZi(wr*) of the product lot* 
of conformal frequency and radius. To leading order at large radius r*, the two independent solutions of equation ((84)) 
are 

sZ™ ~ (a;r*)'* e"*'^''* (propagating for s > 0) , . 
~ (wr*)"" e*"""* (propagating for s < 0) . 
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The modes are labeled in and out because ^'^ oc e *"*sZ™ cx {ujr*y e ''^(*+'' ) corresponds to ingoing modes, while 
(X e~'"'sZ°"* (X (ijjr*)~* e~"^'^*~'' ^ corresponds to outgoing modes. For non-zero spins s, only one of the two 
modes of equations ((85)l is propagating, the one that falls off most slowly at large radius r* . For positive spin +s, the 
propagating mode is the ingoing mode, +s2/" oc {ujr*)^e~'''^^ , while for negative spin — s, the propagating mode is the 
outgoing mode, -sz"^^ oc (ajr*)'*e*"'' . The other of the two modes, for non-zero spin, corresponds to the short-range, 
non-propagating spin-(^s) component of a propagating spin-(±s) wave. 

Since changing the sign of the spin s changes the operator on the left hand side of equation into its complex 
conjugate, it follows that solutions ^z; of opposite spin are complex conjugates of each other. Specifically, ingoing 
solutions are complex conjugates of outgoing solutions of opposite spin: 



(86) 



For vanishing Weyl scalar C = 0, as is being considered in this subsection, the radial operators s'Dy and sl?„ defined 
by equations ([69]) form a consistent ladder of spin raising and lowering operators [cf. the comments in the paragraphs 
containing equations ()65|) and ()75|) ]. As will now be shown, these operators can be used to construct exact analytic 
expressions for the asymptotic solutions gZi by successive raisings or lowerings from the spin zero solutions. 

The spin zero solutions oZi{uir*) of equation ([M)) are proportional to spherical Hankel-Bessel functions h\'''\ujr*) = 
ji{ujr*)±tyi{cor*) = ^Tr/{2ujr*) 1 {Lur*)± iY^+i {cur*)] , 



oz^H^r*) = {iy+^Ljr*h\^>{Ljr*) , qz^'K*) = {-iy+^Lor* h\-> (ujr*) , (87) 

The spin zero solutions Qzl{ujr*), 



which are complex conjugates of each other, in accordance with equation 
equations ([57)) . are normalized so that 



ozf^\u;r*) 



as r 



(88) 



The spin raising and lowering operators gP-u act not on sZi{ujr*) but rather on the full temporal-radial part sZi^i(t, r) 
of = sZujiit, r) sYimiO, 4>), equation dTB]), 



,Z^;(t,r) = {C/n^''^A~^e'''^\zi{Lor*) 



(89) 



■a 

The functions ±sZuji{t,T) are by definition eigenfunctions of the Hermitian operator j-^I?''" j-^Pu, with eigenvalues 
— {l—s){l+s^\)/2 in accordance with equation ([77)) . 



{l-s){l + s + \) 



,Z^z(i,r) = 



The radial spin raising and lowering operators sT^v and sT^u raise and lower the spin s of sZuji{t, r) by one unit 



sVv sZujiit,r) 



{lTs)il±s + l) 



1/2 



iZ^i{t,r) 



(90) 



(91) 



where the normalization factor follows, modulo a phase factor, from the fact that the raising and lowering operators 
gVy and s+i^u are Hermitian conjugates of each other, equation (|70p . and that their product s+i^^u s^v has eigenvalue 
~{l~s){l+s+l)/2, equation 

The reduced radial spin raising and lowering operators sdv and sdu defined by equation (|80p or equivalently (I8ip 
simplify to 



1 



d„ = — ±r* — 



V2 



d 



dr* 



(92) 



Similar to equation ((9T|) . the radial spin raising and lowering operators gdv raise and lower the spin s of sZi{u}r*) by 
one unit 



sdv sZi{ujr*) 



(?Ts)(/±s + 1) 



1/2 



s±l 



zi {ur*) 



(93) 
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There are of course both ingoing solutions sZf^ and outgoing solutions s^™*: the operators sdv raise and lower ingoing 
solutions to other ingoing solutions, and outgoing solutions to other outgoing solutions. The asymptotic behavior 
of the spin zero modes qZi, equation ([55]) . combined with the raising and lowering equations fM]). imply that the 
asymptotic behavior of sZi at large radius is 





.)!] 


[{1 + 






.)! 







1/2 



(94) 



in agreement with the earlier result (|85p . The choice of phase factor — i in equation ([9T|) and consequently (193)) is such 
that ingoing solutions are complex conjugates of outgoing solutions of opposite spin, equation 



The Wronskian of the two solutions 



and 



must be constant, in accordance with the usual rules for 



homogeneous solutions of second order differential equations, and the value of the constant follows from the asymptotic 
behavior (l94t. 



,5 



dr* 



dr* 



2iu} . 



(95) 



IV. PERTURBATIONS TO SELF-SIMILAR BLACK HOLES 



This section sets up the theory of perturbations to spherically symmetric self-similar black hole spacetimes. The 
approach follows the Newman-Penrose formalism [ssl . The section marches successively through perturbations to 
the vicrbein ( §IV Ap . tetrad connections ( §IV Bp . Riemann tensor ( §IV Cp . and Einstein tensor ( piVEp . The results 
will be used in subsequent sections § jVT4 Vill on gravitational, electromagnetic, and scalar waves. 

From the outset, angular perturbations are required to be expandable in spin-weighted spherical hamornics sY, 
and conditions under which this can be done are obtained. Initially, from requiring that all tetrad connections be 
expandable in spin harmonics, it is found that a certain set of gauge conditions is imposed, the spherical gauge, 
equations (|106p . Subsequently, it is found ( ^IVPP that the tetrad connections that impose the spherical gauge 
conditions are precisely those that do not actually appear in any of the wave equations, so that it is consistent to 
expand vierbein perturbations in spin harmonics with no gauge conditions whatsoever. The general rule is established 
that the spin s of any quantity is equal to the sum of the -l-'s and — 's of its covariant components in the Ncwman- 
Penrose formalism. The result is so fundamental that it deserves stating on a line by itself: 

spin ,s ~ sum of +'s and — 's of covariant indices . (96) 

The motivation for working with spin-weighted spherical harmonics, aside from the familiar fact 0, fl8| that the 
angular eigenmodes of wave equations in the Newman-Penrose formalism are spin harmonics, is that spin harmonics, 
or more correctly ,^Yi„, (d, (t>)e~^'^^ , are eigenmodes of the full 3-dimensional rotation group 0(3), parameterized by 
three Euler angles d,(j),x [35|- Spherically symmetric spacetimes are symmetric not merely with respect to rotations 
about the origin, but with respect to the entire 3-dimensional rotation group. Thus it is natural to expect that angular 
perturbations of spherically symmetric spacetimes should be expandable in spin harmonics. 

In this section, a (zero) overscript signifies an unperturbed quantity, while a 1 (one) overscript signifies a 
perturbation. No overscript means the full quantity, including both unperturbed and perturbed parts. An overscript 
is attached only where necessary. Thus if the unperturbed part of a quantity is zero, then no overscript is needed, 
and none is attached. 



A. Perturbed tetrad 



Define the vierbein perturbation ipmn so that the perturbed vierbein is 

e™'' = (C - ^™") en^ , (97) 

with corresponding inverse 

e'"M = (C + ¥'n")e'V . (98) 
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Since the vierbein perturbation ipmn is already of linear order, its indices can be raised and lowered with the 
unperturbed metric, and transformed between tetrad and coordinate frames with the unperturbed vierbein. The 
perturbation ipmn can thus be regarded as a tensor field defined on the unperturbed background. The perturbation 
Lpmn has 16 degrees of freedom, but only 6 of these degrees of freedom correspond to real physical perturbations, 
since 6 degrees of freedom are associated with arbitrary infinitesimal changes in the choice of tetrad, which is to 
say arbitrary infinitesimal Lorentz transformations, and a further 4 degrees of freedom are associated with arbitrary 
infinitesimal changes in the coordinates. 

Under an infinitesimal tetrad transformation, the covariant vierbein perturbations ^Pmn transform as 

(99) 

where emn is the generator of a Lorentz transformation, which is to say an arbitrary antisymmetric tensor. Thus the 
antisymmetric part ifmn — y^nm of the covariant perturbation ipmn is arbitrarily adjustable through an infinitesimal 
tetrad transformation, while the symmetric part ipmn + y^nm is tetrad gauge-invariant. It is easy to see when a 
perturbation to a quantity is tetrad gauge-invariant: it is tetrad gauge- invariant if and only if it depends only on the 
symmetric part of the vierbein perturbation, not on the antisymmetric part. 

Under an infinitesimal coordinate gauge transformation, the coordinates change by an infinitesimal shift Ae^ 
(the factor A is incorporated so that is dimensionless) 

a;'" -> -f Aef" . (100) 

Because the dimensionless shift e'* is, like the vierbein perturbations (fmni already of linear order, its indices can be 
raised and lowered with the unperturbed metric, and transformed between coordinate and tetrad frames with the 
unperturbed vierbein. Thus the dimensionless shift can be regarded as a vector field defined on the unperturbed 
background. Under an infinitesimal coordinate gauge transformation (jlOOp . the vierbein perturbations ipmn transform 
as 

^mn '^mn — {inkdm + ^nkm — ^nmk) Ae*' , (101) 

in which = e'^^ e'^ are the (dimensionless) tetrad components of the coordinate shift. Explicitly, in the Newman- 
Penrose tetrad the perturbations (pmn transform under infinitesimal coordinate gauge transformations (jlOOp as 
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Perturbations can be classified into polar modes, which do not change sign, and axial modes, which change sign, 
when the azimuthal angular tetrad axis 70 is flipped in sign, equation ([M]) . In the Newman-Penrose formalism, 
flipping the azimuthal axis 7^ is equivalent to swapping covariant angular indices -!-<->—, equations ()55|) . which, for 
quantities that are real in an orthonormal tetrad, is in turn equivalent to taking the complex conjugate. Thus in the 
Newman-Penrose formalism, polar modes effectively correspond to real perturbations, while axial modes correspond to 
imaginary perturbations. It is a fundamental feature of the Newman-Penrose formalism not to separate out the polar 
(real) and axial (imaginary) perturbations, but rather to treat them together as a combined complex perturbation. 
The advantage of the Newman-Penrose approach is that the radial part of the wave operator in the wave equation for 
non-zero spin components of the Weyl tensor proves to be the same for both polar and axial modes, revealing their 
fundamental unity [see for example the discussion in the paragraph containing equation (jl78p ]. 

It should be commented that the general strategy for linearization, which has already been invoked in obtaining 
equations (|102p from equations (llOip . and which is used throughout the rest of this paper, is to discard products of 
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quantities that are already first order, since the result is of second order, and to replace directed derivatives of first 
order quantities by the spin raising and lowering operators gT^rm equations (|6ip and (j69p . as follows: 



dv +sTuvv + ar_^+ + 6r+„_ ^ A^+'-'^-^s^v ^-2-^+'^+'' (i03a) 

u vuu +ti- -tl+ « 

9±-sr_++ A^^isPi . (103b) 



+- 



B. Perturbed tetrad connections 



This subsection presents expressions for the perturbed tetrad-frame connections Tkmn in terms of the vierbein 
perturbations (fmn , equations (|104p . Much of the subsection is concerned with showing that a certain set of gauge 
conditions, spherical gauge, equations (|106p . emerges if it is required that the relations (|104p between the tetrad 
connections and vierbein perturbations become algebraic upon resolution into spin- weighted harmonics, that is to 
say, if each (perturbation to a) tetrad connection and each vierbein perturbation is required to be individually 
expandable in spin-harmonics of some definite spin. The next subsection, ijIV C[ continues the theme, finding that 
spherical gauge emerges a second time if it is required that the relations between the Riemann tensor and the tetrad 
connections become algebraic on resolution into spin- weighted harmonics. Ultimately, in §IVD| the need for spherical 
gauge is traced to just 4 offending tetrad connections [the ones involving boxed terms in equations (|104p below], and 
it is pointed out that these 4 connections are precisely the ones that do not actually appear in the wave equations in 
the Newman-Penrose formalism. Thus not only are the spherical gauge conditions not required, but in fact the wave 
equations in the Newman-Pcnrose formalism can be resolved into spin-weighted spherical harmonics with no gauge 
conditions at all. 

Three important lessons emerge from the analysis of this and the next two subsections. First, it is consistent to 
expand vierbein perturbatioris iOmn in spin- weighted harmonics from the outset: there is no need for the more general 
system of tensor harmonics 28, Sa. [gsI [93|. Second, the wave equations in the Newman-Penrose formalism become 
algebraic upon resolution into spin-weighted harmonics without any gauge conditions; that is, one is free to choose any 
gauge whatsoever. However, if a gauge is called for, then the spherical gauge conditions (jl06p are, though not required, 
nevertheless natural, since it is only in spherical gauge that each and every tetrad connection can be expanded in spin 
harmonics of definite spin. Subsequently, W C[ it will be found that spherical gauge proves to provide a particularly 
elegant gauge for expressing the spin-±2 wave equations that govern the propagating components of gravitational 
waves. Third, the spin-weight of any object is given by the sum of the -l-'s and — 's of its covariant angular indices in 
the Newman-Penrose formalism, as already asserted by rule ([M)) . 

To linear order in the vierbein perturbations if mm the perturbed tetrad- frame connections Tkmm computed from the 
fundamental formula psp . are (the boxed terms in the following equations are terms that break algebraic resolution 
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into spin-weighted harmonics, as discussed below) 
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In the unperturbed self-similar background, the non-vanishing Newman-Penrose components of the tetrad-frame 
connections Tkmn reduce to 
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The tetrad connections Tkmn are, by construction, all coordinate gauge- invariant, because they are defined, 
equation (|77|) . purely in terms of tetrad frame quantities, independent of the choice of coordinates. However, for tetrad 
connections whose unperturbed part is finite (non-zero), the split between the unperturbed and perturbed parts of 
Tkmn is not coordinate gauge- invariant. Moreover, the tetrad connections Tkmn are not tetrad gauge- invariant. 

It is a striking fact that equations (|104p for the tetrad connections Tkmn become algebraic upon resolution into 
spin-weighted harmonics each perturbation tpmn and each connection Tkmn being expanded in harmonics of 
appropriate spin s, but only provided that certain gauge conditions are imposed, namely that the perturbations ipmn 
vanish whenever the right index n is an angular index, that is, n = -I- or — : 



1^++ = = 



(106) 



These 8 gauge conditions can be termed "spherical gauge" , because physically they correspond to choosing a set of 
perturbed coordinate and tetrad frames with respect to which the spherical symmetry of the background is preserved 
to the greatest extent. As argued immediately below, spherical gauge can be accomplished by 3 coordinate and 5 
tetrad gauge transformations. 

To see how the spherical gauge arises mathematically from equations (|104p , consider for example the perturbation 

<fv+, which appears operated on by the angular operator +iD_ in equation (|104bp for ^ ( r_t,+ -|-r+,;_ j , and by ^iDj_ in 

\ +11- -M+/ 

equation (|104kp (the boxed term containing (pv+) for r+„+, but by ^iD_ in equation (jl04dp for r_+„ — i( r_„+— r+t,_ ) . 

-u- +-U \ +u- -u+J 

Whereas and ^iD^ expect to operate on a spin-harmonic sY of spin-weight s = -f 1 (in accordance with the 
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prefix s on ^2?^), by contrast _iD_ expects to operate on a spin-harmonic of spin- weight s = — 1. This discrepancy in 
spin- weights can be accomodated only if vanishes. Similar discrepancies occur between equations (|104bp or (jl04k[) 
and (jl04d|l for (fv-, 'Pu+, and (fiu^, and again these can be accomodated only if all these perturbations vanish. The 4 
perturbations fv+, fv-, fu+i and <^«_, can be made to vanish by 4 infinitesimal tetrad gauge transformations. It is 

consistent to point the finger of blame at the boxed terms in equation (Il04dp for r_+^ — ^ I r_^+ — T+v~ ] ■ The general 

V +u- -u+J 

rule that emerges below, equations piOp and (jllip . is that the spin- weight s of a quantity equals the sum of the 
+'s and — 's of its covariant indices, a rule already anticipated in equation ([96|) . The boxed terms in equation (|104dp 
violate the rule. 

Similar discrepancies occur for the perturbations ip++, (p--, </5+-, and ip^^. For example, ip^^ appears operated on 
by in equation (|104ep (the boxed term containing </?++) for r_++, which requires that ipj^.+ be a harmonic of 

spin-weight s = +1, whereas equation (|104k[l for T^y^ [or alternatively equation (jl041| for r+u+J requires that ip+^ be 
a harmonic of spin-weight s — +2. Consistency requires (^++, and similarly iy9__, to vanish. In the case of (/3+„ and 

(p-+, the cot term multiplying these perturbations in equation (|104ep (the boxed term containing ip+-) for r_++ 

. . . ~+ 

destroys algebraicity of the spherical transform unless and (p_^ vanish. Note that it is fine that the unperturbed 

value of r_++ be proportional to cotfl; it is the perturbations, not the unperturbed values, that are to be expanded in 

H 

spin- weighted spherical harmonics. Under an infinitesimal coordinate gauge transformation, ip_f_^ and ip__ transform 
according to equation (jl02hp . and can thus be set to zero by a suitable choice of and e'^ . Then ip^_ + iy9_ + , which 
transforms, in accordance with equation ( |102gp , as 

<y5+- + <P-+ + <^-+ - (2 Pve" + 2 /3„e" + +iV_ + e+) , (107) 

can be set to zero by a suitable choice of /3 • e = /3^e^ + /3«e". Finally, — can be set to zero by an infinitesimal 
tetrad gauge transformation. 

The conformal factor A transforms under an infinitesimal coordinate gauge transformation as 

A^A + Ae"'d,nA = A{l + f3-e). (108) 

Thus a gauge choice oi (3 ■ e, as accomplished by the vanishing of (^+_ + + , equation (|107p . is equivalent to fixing 
the gauge of the conformal factor A. 

Monopole and dipolc harmonics, sYim with / = or 1, do not require special treatment. For example, equation (jl02hp 
constitutes an algebraic relation between the Zm'th coefficient of the expansion of (^++ in spin ±2 harmonics ±2^m, 

and the Zm'th coefficient of the expansion of e'f in spin ±1 harmonics ±iYim- For the dipole harmonic, I — 1, all 
terms in equation (jl02hp are identically zero, since there is no spin ±2 dipole harmonic, while yields zero when 

acting on any spin ±1 dipole harmonic ±iYim. Thus for the dipole harmonic, equation (Il02hp docs not represent a 
coordinate gauge freedom, but rather it vanishes identically. 

In spherical gauge, two infinitesimal gauge freedoms remain, 1 coordinate freedom and 1 tetrad freedom. The 1 
coordinate freedom corresponds to an infinitesimal adjustment of the conformal radius r [the gauge of the conformal 
factor A being fixed in spherical gauge, in accordance with the comments following equation pOSp ], while the 1 
tetrad freedom corresponds to an infinitesimal Lorentz boost in the radial direction, that is, to an adjustment of the 
antisymmetric part of (p^u- 
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In spherical gauge, equations (I106p . the tetrad-frame connections Tkmn, equations (|104p . reduce to 
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As intended, these equations become algebraic upon spherical transform into spin-weighted harmonics, with the 

1 

perturbations v?„m and the perturbed connections Ti^mn being proportional to spin-harmonics sY as follows: 
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It is apparent from equations (jllOp and (jllip that the spin-weight s of any term is just equal to the sum of the +'s 
and — 's of its covariant indices, in agreement with the previously asserted general rule (j96p . 

It is useful to record that in spherical gauge the directed derivatives of the conformal factor A are, to linear order. 



dvA = ^A{r-v+ + T+V- 

u \ +u- -u+ 
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C. Perturbed Riemann tensor 



Much of this subsection is concerned with demonstrating that the spherical gauge, equations (jl06p . emerges a 
second time, independent of the arguments of the previous subsection, i 3IVB( if it is required that the linearized 
equations (|113p - (|117p for the tetrad components Rkimn of the Riemann tensor in terms of the tetrad connections 
rfcrnra bccomc algcbraic upon resolution into spin- weighted spherical harmonics. Essential to the argument is that 
every tetrad connection be individually expandable in spin harmonics. In the next subsection, jjlVDl it will be shown 
that indeed it is the condition on the tetrad connections that imposes spherical gauge: without the condition on the 
connections, no gauge conditions are imposed. 

General expressions for the Newman-Penrose components of the Riemann tensor Rkimn in terms of the tetrad 
connections Tkmn and their directed derivatives are given in the notation of the present paper in the Appendix, 
equations (jAip . To linear order, but in a general gauge, not in spherical gauge, the Newman-Penrose components of 
the Riemann tensor, equations (jAip . reduce in accordance with the linearization procedure described in the paragraph 
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containing equations (|103p to the expressions (|113P " (|117p given below. For clarity, the equations are subdivided into 
batches characterized by the dependence of their perturbations on spin s and on infinitesimal tetrad gauge (Lorentz) 
transformations. As will be demonstrated below, the spin s of any component Rumn of the Riemann tensor is equal 
to the sum of the +'s and — 's of its covariant indices, a statement that will be found in jjlVDI to be true in any gauge, 
not just spherical gauge. In the spherically symmetric unperturbed background, all components of zero spin except 
Rvu+~ are are in general non- vanishing, while components of non-zero spin necessarily vanish. 

The boxed terms in the expressions (|113p ~ (|117p are either (a) non-conforming terms that break algebraic resolution 
into spin harmonics and hence lead to spherical gauge, as discussed below, or (b) terms proportional to r_+„ or r+_„, 
expressions for which in terms of vierbein potentials were previously found, equation (jl04dp . to contain non-conforming 
terms. 

Perturbations of the following 5 components of the Riemann tensor are proportional to spin-0 harmonics qY , and 
are tetrad gauge-invariant with respect to all 6 arbitrary infinitesimal tetrad transformations. Because of the Jacobi 
identity (|114p . the 5 components represent only 4 independent degrees of freedom. The components Gyu and 
of the Einstein tensor, §IVE| and the spin-0 component Co of the complexified Weyl tensor, W A\ depend on these 
Riemann components: 
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The last three Riemann components of equations (jll3p are related by the only non-trivial occurrence of the Jacobi 
identity: 
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Perturbations of the following 2 components of the Riemann tensor are proportional to spin-0 harmonics qY, 
but are not tetrad gauge-invariant, since they vary under an infinitesimal tetrad transformation associated with the 
antisymmetric part of fvu, that is, under an infinitesimal Lorentz boost in the radial direction. The components Gm 
and Guu of the Einstein tensor are proportional to these Riemann components: 
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The following 8 components of the Riemann tensor are proportional to spin-±l harmonics ±iY . They are not tetrad 
gauge- invariant: components proportional to _^iY vary under transformation of the antisymmetric parts of Lp^^ and 

while components proportional to _iY vary under transformation of the antisymmetric parts of Lpy_ and 
They are however tetrad gauge- invariant with respect to an infinitesimal radial Lorentz boost, the one tetrad gauge 
freedom remaining in spherical gauge. The components Gt,+ , Gv_, Gu+, and G„_ of the Einstein tensor, and the 
spin-±l components C±i components of the complexified Weyl tensor, depend on these Riemann components: 
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(116a) 
(116b) 

(116c) 
(116d) 
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The following 6 components of the Riemann tensor are proportional to spin-±2 harmonics ±2Y. They are tetrad 
gauge-invariant with respect to all 6 arbitrary infinitesimal tetrad transformations. The components G^^ and G__ 
of the Einstein tensor, and the spin-±2 components C±2 of the complexified Weyl tensor, depend on these Riemann 
components: 

= -A-~^±iVur+^++A-\iV±T+vu, (117a) 

V —U— -UV 

Rv+v+ — — ±iDv A j^v+ ^ A ^^iD±Tj^vv , (117b) 

u—u— u —u— -uu 

Ru+u+ = - ±iVuA-^T+u++ A-^^{D^T+uu . (117c) 

V — V— V —V— —vv 

In the unperturbed background, the non-vanishing components of the tetrad frame Riemann tensor are 

2 2 

{DuA^^K , (118a) 

R+_^_ = -2A-^M , (118b) 

Rv+u- = Rv-u+ - - A-^^iU, A-^Pu = - A-\iVu A-'f3y , (118c) 

= -A±iVvA-^Pv. (118d) 

U-U+ u u 

In terms of the dimensionless quantities R, P, P±, and F defined by equations (l39|) and related to the proper density, 
radial pressure, transverse pressure, and energy flux through Einstein's equations ([42]) . the unperturbed Riemann 
components are 

Rvuvu = (R- P + 2P±-2M)/A\ (119a) 
= -2M/A2 , (119b) 







Rv+u- = Rv-u^ = -^{R-P- 2M)/A^ , (119c) 
Rv+v^ = UR + PT'2F)/A^ . (119d) 



U — U + 



It is remarkable that equations (I113p - ()117p for the perturbed Riemann tensor in terms of the tetrad connections 
become algebraic upon resolution into spin-weighted harmonics, each connection Tkmn and each Riemann component 
Rkimn being expanded in harmonics of appropriate spin s, but only in a particular gauge, and this is the same spherical 
gauge, equations (|106p . as obtained in WV Bl from the requirement that the expressions (|104p for the connections in 
terms of the vierbein perturbations become algebraic upon resolution into spin- weighted harmonics. 

The argument that equations (|113p - (|117p imply spherical gauge runs as follows. First, if a connection Tk^n appears 
operated on by a spin raising or lowering angular operator sTJ± anywhere in equations (|113p - (|117p . then it follows 
immediately that the connection Tkmn must have spin s. Inspection of equations (|113p - (|117p shows that all 16 of the 
connections that vanish in the unperturbed background do indeed appear somewhere acted on by an angular operator 
si'i, and the spin s of those 16 connections then follows immediately. As for the 8 connections whose unperturbed 
values are finite (non-zero) , all 8 of those connections appear acted on by an angular directed derivative d± somewhere 
in equations (|113p - (|117p . In the case of the 6 connections whose unperturbed part is a conformal factor times a purely 
radial function, that is, r^„„, Tuw, T-d+j r+^u-: r+u-: and r_„_|_ [every case except r__|_+ and r_|___, which will be 

considered below, equations (jl22p ]. these terms can be split into parts depending on the unperturbed connection F 

1 

and its perturbation F as follows [an angular directed derivative d± acting on a first order quantity can be replaced 
by A~^s'^± in accordance with the general rule (jl03bp . and in every case in equations pi3p - (|117p where d± acts on 
a connection whose unperturbed value is finite, the spin s is zero]: 

5± F = (ifi+u dv + ip+v du)T + A-^oV^ V . (120) 

1 

Here it is necessary to assume not only that the perturbation F has a definite spin s, but also that each vierbein 
potential ^^v, ¥'-i>, and iy9_|_„ individually has definite spin s. One way or another, all 24 tetrad connections 

appear acted on by an angular operator somewhere in equations (|113p - (lll7p . and the net consequence is that the spin 
s of every connection [with the possible exception of F_++ and F+__, considered explicitly below, equations p22p ]. 
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or of its perturbation if the unperturbed value is finite, is equal to the sum of the +'s and — 's of its covariant 
indices, in agreement with the earlier conclusion, equations (jllip . The spin s of each of the Riemann components 
Rkimn then follows, and again the result is that the spin s of the Riemann component, or of its perturbation if the 
unperturbed value is finite, is equal to the sum of the +'s and — 's of its covariant indices, in agreement with the 
general rule (|96p . Radial directed derivatives dv acting on connections whose unperturbed values are finite also appear 
variously in equations (|113p - (|117[) . For the 6 connections whose unperturbed part is a conformal factor times a purely 
radial function, the perturbations to the radial directed derivatives of these can be split into parts depending on the 

1 

unperturbed connection V and its perturbation F as follows: 



du)T + dvT 

II. ^ lil) II. 11.11. l> ' II. 



(121) 



Here again it is necessary to assume that the vierbein perturbations ^}uv^ Vw, and (puu each have definite spins. In 
combination with the constraints on i^+i,, ^P-u, V^vi and one concludes that the spin of each vierbein perturbation 
</'mn whose right index is u or w is equal to the sum of the +'s and — 's of its covariant indices, again in agreement 
with equations (|110p and in accordance with the general rule . 

So far it has been concluded that, if it is assumed that (the perturbation of) every tetrad connection Tkryim and 
every tetrad Riemann component Rkimm and in addition every vierbein potential (pmn whose right index n is u or u, 
is expandable in spin harmonics sY of definite spin s, then equations (|113p - (|117p for the Riemann tensor in terms 
of the tetrad connections require that the spin of each object is equal to the sum of the +'s and — 's of its covariant 
indices. 

The next step in establishing spherical gauge is to identify in equations pi3p - (|117p non-conforming terms, those 
that do not have the correct spin. These terms call attention to themselves by involving F_++ or F+__, because the 
latter connections are proportional to cot in the unperturbed background, and misfitting factors of cot destroy 
algebraicity of the spherical transform into spin-weighted harmonics. That is, cot 9 times a spin-weighted spherical 
harmonic sYim is not another pure spin-weighted spherical harmonic. Thus the boxed terms (F+„„ — F+„t,)F_|___ and 

both be zero. Similarly the boxed 



\^ — uv 

terms F+ 



F__|_+ in equation (|113ep require that F 



F+„„ and F_ 



F+t,+ in equation (|116cp . and F+__ F+m+ in equation (jll6dp . require that F^ 



and F+tj+ all 



be zero. Other non-conforming terms in equations (|113p - (|117l) . boxed for clarity, are those where F^^^, or F+_„ is 
acted on by a directed derivative, instances of which occur in equations (I113bp . (|116cp . and (|116dp . These terms are 
explicitly 



(a±+F 
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^ (- 1 + <^++ + <^+-) + ^ (o2?± - o2?T + A-^^i'Di r+- , (122a) 



-+ 
cot ( 



(A^oVvA-^(p+--±iV±ipv-) +QVvA~^r^++, (122b) 

II. I II \ II. I ^ ' 

. (122c) 



1 cot 1 1 
-((pu+ + ipu-) H ^(^%X>uA"V+- - ±i'D±ipu-) +o'DuA-^r 

2 V- V+' y'2 ^ 1- U+'J V 



Each of the terms multiplying cot0 in equations (|122p destroys algebraicity of the spherical transform into spin- 
weighted spherical harmonics, and from these alone all of the spherical gauge conditions (|106p necessarily follow. In 
spherical gauge, T^yu — F+k„, T_uv — F_^„, F+^+, F_„_, F_i,_, and F+„+ all vanish as required, and equations (|122p 
reduce to 



(dv 

^ U 
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H — u 
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+u 
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u+ I 



(123a) 
(123b) 
(123c) 



which now conform to being algebraic upon transform into spin-weighted harmonics, provided that the spins of the 

1 1 ___ 

perturbations F_+^. and F^.__ are -1-1 and —1 respectively, in agreement with the general rule (j96p that the spin of 

an object equals the sum of the +'s and — 's of its covariant indices. 

This completes the demonstration that, as claimed, equations (|113p - (|117p for the Riemann tensor in terms of the 

tetrad connections become algebraic upon transformation into spin- weighted harmonics in and only in spherical gauge, 

equations p06p . 
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D. To gauge or not to gauge? 



Much of the previous two subsections, iji jlVBI and llV C[ was concerned with showing that a special gauge, spherical 
gauge, equation p06p . emerges from demanding either that every vierbein perturbation (fmn and every tetrad 
connection Tkrnn be individually expandable in spin- weighted harmonics, or alternatively that every tetrad connection 
Tfemn and every tetrad Riemann component Rkimn be individually expandable in spin-weighted harmonics. 

But what if instead one merely requires that each vierbein perturbation fmn and tetrad Riemann component Rumn 
(but not necessarily each tetrad connection) be individually expandable in spin-weighted harmonics? As will now be 
shown, this requirement is satisfied with no gauge conditions at all. 

Inspection of equations (|104p for the tetrad connections in terms of the vierbein perturbations shows that only 4 
of the connections contain non-conforming terms (boxed terms) that break algebraic resolution into spin harmonics, 
namely r_+u, equation (jl04dp . and r_++, equation (|104ep . In effect, it is the requirement that these 4 connections be 

expandable in spin harmonics that imposes the spherical gauge conditions (jl06p . The offending 4 tetrad connections 
appear in just 6 components of the Riemann tensor, the boxed terms in the expressions for Ryu+-, R+-v+, 

_ -+u- 

and R+-U+, equations (I113bp . pi3ep . (|116cp . and (IllOdP . However, 5 of these 6 Riemann components have alternative 

-+v- 

expressions in terms of tetrad connections and their derivatives, and in each case the alternative expression contains 
only conforming terms (no boxed terms). Since the alternative expressions in each case must yield the same expression 
in terms of vierbein perturbations, it follows that the more fundamental expressions for the 5 Riemann components 
in terms of vierbein perturbations contain only conforming terms. This leaves only 1 potentially non-conforming 
Riemann component, R^_^_. The non-conforming (boxed) terms in this component, equation (|113bp . are expanded 
in equation p22ap . but the two parts combine to give 

(a+ + r_++)r+__ + (9_ + r+__)r_++ = a-^ [- 1 + + + ^iV- ^v^{^^_ + ^_+) 

+22?- <P++ - -22?+ (124) 

all terms of which now conform. [For later reference, the left hand side of equation (|124p equals 2A^'^M where M is 
defined by equation p90p ]. 

Thus no gauge condition is imposed if it is required only that each vierbein connection (pmn and each tetrad 
Riemann component Rkimn be expandable in spin harmonics (without any conditions on the tetrad connections). 

The statement of the previous paragraph can be put another way. If the vierbein perturbations fmn are expanded 
in spin- weighted harmonics with spin given by the general rule (1961) , then it automatically follows that each Riemann 
component Rkimn will be a sum of spin- weighted harmonics with spin given by the general rule (|96p . 

Since it is now evident that the spherical gauge conditions (|106p follow precisely from requiring that the tetrad 
connections be expandable in spin harmonics, the question arises as to whether or not it is necessary to impose 
these conditions in order to carry through the wave equations in the Newman-Penrose formalism. Recall that 
only 4 of the tetrad connections were non-conforming in the first place, namely F-^^, equation (I104dp . and F-+4., 

equation (|104ep . Cunningly, it turns out that the linearized equations governing gravitational, electromagnetic, and 
scalar waves conspire to involve all of the tetrad connections except the 4 non-conforming tetrad connections! Notably, 
the expressions (|157p for the Weyl currents Jimn involve all but the 4 non-conforming tetrad connections. The 4 non- 
conforming tetrad connections seemingly appear in the expressions (I154p . (jl62p . (|239p . and (|245p for the differential 
operators s\n, sA'„, sSm, and but everywhere that these operators actually occur in the relevant linearized 

equations, namely equations (|156p . (|165p . (|24ip . (|248p . and (|286p . the spin s of the differential operator is such that 
the contribution of non-conforming tetrad connections vanishes. 

Thus it can be concluded that the wave equations in the Newman-Penrose formalism can be expanded in spin- 
weighted harmonics, and that this can be done without any gauge conditions. Put another way, one is free to choose 
any gauge that is convenient. 
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E. Perturbed Einstein tensor 



In a Newman-Penrose tetrad, the Einstein tensor Gmn is related to the Riemann tensor Rkimn by 

(125a) 
(125b) 
(125c) 

(125d) 

(125e) 

(125f) 

Because of the variety of expressions for the individual Riemann components, equations ()113|) - (I117|) . there are several 
alternative ways to express each Einstein component in terms of tetrad connections and their derivatives (though only 
one way to write each of G++ and G__), and there is no pressing reason to write out these expressions. 

In the unperturbed self-similar background, the non- vanishing Newman-Penrose components of the Einstein tensor 
Gmn are, in terms of the dimcnsionlcss quantities R, P, P±, and F defined by equations 
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G^u = {R- P)/A^ , (126a) 

G+_ = 2P_i/A2 , (126b) 

G™ = {R + P- 2F)/A^ , (126c) 

Gun = {R + P + 2F)/A^ . (126d) 
It is convenient also to define an auxiliary dimensionless unperturbed quantity E by 

^= iA2(G,„ + G+_) = i(i?-P + 2P^) , (127) 

since this quantity occurs several times in the wave equations, notably equations p72p and (|199p . 

It follows from equations (I125P that the perturbation of each component of the tetrad-frame Einstein tensor follows 
the Riemann tensor in conforming to the general rule ([M)) that its spin- weight s equals the sum of the -l-'s and — 's of 
its covariant indices: 

Gvu oc Gvv OC 0^ , (128a) 

-f- uu 

Gv+ (X G„+ cx ±ir , (128b) 

U— V- 

G++ (X ±2Y . (128c) 

Being tetrad- frame quantities, all the components Gmn of the tetrad- frame Einstein tensor are automatically 
coordinate gauge- invariant, although the split between the unperturbed and perturbed parts of components whose 
unperturbed parts are non-zero, namely Gyu, G+_, G^, and Gum is not coordinate gauge-invariant. 

The components Guv, G+_, G++, and G__ are tetrad gauge-invariant with respect to all 6 arbitrary tetrad (Lorentz) 
transformations. The remaining components are not tetrad gauge-invariant. The components Gyy and G„„ vary 
under transformation of the antisymmetric part of ifyu-, that is, under an infinitesimal radial Lorentz boost, the one 
tetrad gauge freedom left in spherical gauge. The components G„+ and G„+ each vary under transformation of the 
antisymmetric parts of ipu^ and fu+, while G^^ and Gu- each vary under transformation of the antisymmetric parts 
of (fiv- and ipu-- The following quantities are tetrad gauge-invariant, but not coordinate gauge- invariant: 



Gvv(l + ifiuv — (fivu) , (129a) 

uu ^ vu uv^ 

Gv+ + - V+v)iGvu + G+_) + h(ipu+ ~ (p+u)Gvv , (129b) 

u- ^ u- -u' V- -v' uu 

Gv- + ^{ifiv- - ip-v)iGuu + G+_) + ^{fu- - ip-u)Gvv . (129c) 



Note that the product G^Guu is both coordinate and tetrad gauge-invariant. 

There are just two components of the Einstein tensor, the spin-±2 components G++ and G_„, that are both 
coordinate and tetrad gauge-invariant, and also vanishing in the background, and therefore constitute physical 
perturbations that are unchanged by any gauge transformation. 



28 



In the limiting case of the Reissner-Nordstrom geometry, the only components of the Einstein tensor that are 
non- vanishing in the unperturbed background are G„„ = Gj^_. In this case the components Gy^ and Guu are also 
coordinate and tetrad gauge-invariant and vanishing in the background. 

Finally, in the case of vanishing background energy-momentum tensor, the Schwarzschild geometry, all perturbations 
of the Einstein tensor are coordinate and tetrad gauge- invariant. 



V. GRAVITATIONAL WAVES 



This section derives the equations that describe gravitational waves in perturbed self-similar black hole spacetimes. 
Subsection lV Al characterizes perturbations of the Weyl tensor, WBI presents the equations that govern the evolution of 
those perturbations, and subsequent subsections apply those equations to derive gravitational wave equations, starting 
with W C[ which presents the wave equations governing the propagating (spin-±2) components of gravitational waves. 
Subsection IV Dl discusses how the Weyl evolution equations relate different spin components of gravitational waves to 
each other, so that the different spin components fluctuate in harmony rather than independently. Subsections IV El 
and IV Fl focus on the special cases of monopole and dipole modes. In the final two subsections IV Gl and IV Hi of this 
section, the spin-0 wave equation in the Newman-Penrose formalism is used to derive axial and polar spin-0 wave 
equations. 



A. Perturbed Weyl tensor 



Gravitational waves are most naturally described by the Weyl tensor Ckimn, the totally trace- free, or tidal, part of 
the Riemann tensor Rkimm equation psp . Like the Riemann tensor, the Weyl tensor is a symmetric bivector matrix, 
that is, a symmetric matrix of antisymmetric tensors. 

As is familiar from the geometric algebra [26| . a bivector is a 6-component object which, in an orthonormal frame, 
has a natural complex structure (indeed, a natural complex quaternionic structure; but the quaternionic structure 
goes beyond what is needed for the present paper). The real, or "electric" (E), part of the bivector changes sign 
under spatial inversion, while the imaginary, or "magnetic" (B), part does not change sign under spatial inversion. 
In the geometric algebra, it is natural to treat a bivector as a single complex object E + iB. 

In an orthonormal tetrad, the Weyl tensor, being a symmetric matrix of bivectors, can be organized as a 2 x 2 
matrix of 3 x 3 blocks, with the structure 



EE EB 
EB BB 



(130) 



where EE, EB, and BB signify 3x3 matrices with components of the indicated type, E for electric, and B for 
magnetic. In view of the bivector structure, it is natural to define a complexified version Ckimn of the Weyl tensor by 



Cklrnn ^ - + -S.l^A (S^j:, + -e,r.n'-') Cp.rs , (131) 



2' 

where eupq is the totally antisymmetric tensor, here normalized so that etrs^ = — 1. The overall factor of 1/4 on the 
right hand side of equation (|13ip is introduced because then the operator P^^ = i + ^Ski^"^) is a projection 
operator, satisfying — P. The complexified Weyl tensor Ckimn is, like the Weyl tensor, totally traceless. If the 
Weyl tensor Gkimn is organized according to the structure (|130p . then the complexified Weyl tensor Ckimn has the 
corresponding structure 

i Q ^^yEE-BB + 2^EB) . (132) 

In tensor language, the complexified Weyl tensor satisfies the symmetries 

Gki "Ckirs ■ (133) 

Thus the independent components of the complexified Weyl tensor Ckimn constitute a 3 x 3 complex symmetric 
traceless matrix EE—BB+2iEB, with 5 complex degrees of freedom. 

The complexified Weyl tensor Ckimn remains unchanged under spatial inversion, because spatial inversion not only 
transforms E — > —E and B ^ B, but also changes the sign of the antisymmetric tensor Skipq, thereby effectively 
changing the sign of i in EE—BB+2iEB. 
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In a Newman- Penrose tetrad, the non- vanishing components of the complexified Weyl tensor Ckimn, and their 
expressions in terms of the Weyl and Riemann tensors, are 

Cq — Cuvuv — ^_ Cuv-\-- — C^^^u — 2 {^uvuv ^ — 2 (^H h- ^~ Cuv + ~^ — Cy^_u 

— ^ {Ruvuv + — h- + 2 Ry^_u + 2 Ruv+-) , (134a) 

C±l = — (^-+1)+ = CvUV+ — C-+tf+ — h (RvUV+ + ; (134b) 

UVU- +-U- UVU— +—U— uvu- +~u- 

C±2 = Cv+v+ — Cv+v+ — Rv+v+ ■ (134c) 
u — u~ u — u— u — u~ 

Notice that the only non- vanishing components of the complexified Weyl tensor Ckimn are those with bivector indices 
uv, H — , v+, or Components Ckimn with bivector indices v— or u+ all vanish: 

Cv~^u — Cvuv- — C+-V- — Cv-v- — , (135a) 

UVU+ -+U+ u-fu + 

Cv—v+ — Cv-u- — Cv-v+ — Cv-u- = , (135b) 

U+U- U+V+ U+U- U+V+ 

in which the top and bottom lines, equations (jl35ap and (jl35bp . collect components Ckimn for which the parent Weyl 
tensor Ckimn respectively does not, and does, also vanish. 

In the unperturbed self-similar background, only the spin-0 component Co of the complexified Weyl tensor is 
non-zero, and it equals 



Co = C/A^ (136) 



where C is the unperturbed Weyl scalar given by equation p4|) . 

The Cs defined by equations (I134p can be referred to as the spin-s components of the complexified Weyl tensor, 

since Cg (or in the case s = the perturbation Co) is proportional to a harmonic of spin- weight s 

1 

Cs cx sY (137) 

as follows from equations p34p and the fact, demonstrated in mV C\ that each component of the Riemann tensor 
conforms to the general rule ([96]) that its spin equals the sum of the -l-'s and — 's of its covariant indices. As expounded 
in WCI below, the spin-±2 components C±2 of the complexified Weyl tensor describe propagating components of 
gravitational waves: C+2 describes the propagating component of ingoing gravitational waves, while C_2 describes 
the propagating component of outgoing gravitational waves. In more traditional notation [H, [63, [H, [8§| , Cg = 

so that C^2 is ^'o, while C_2 is ^'4 (Chandrasekhar's [l§| metric signature H is opposite to the present paper's 

— \ — \-+, so his Riemann tensor is opposite in sign to the present paper's). 

The complexified Weyl tensor Ckimn, even though defined, equation ()13ip . as a projection of the Weyl tensor Ckimn, 
nevertheless retains all the degrees of freedom of the Weyl tensor. This can be demonstrated by resolving Cg into 
polar (p) and axial (a) parts 

Cg = C^P^ + Ci'') , (138) 

defined by the property that polar and axial parts respectively do not and do change sign when the azimuthal angular 
tetrad axis 7^ is fiipped in sign, equation The angular conjugate of Cg is 

C* = Ci"^ - Ci°) . (139) 

As remarked after equation (|55p . the effect of angular conjugation is to swap angular indices -|- <-!■ — on tensors (such 
as the Riemann tensor) and other natively real objects in a Newman- Penrose tetrad (such as the tetrad connections); 
but angular conjugation does not swap angular indices on a complexified object such as the complexified Weyl tensor 
Ckimn- Thus the angular conjugates C* of the spin components of the complexified Weyl tensor are given in terms of 
the Weyl and Riemann tensors by the angular conjugates of equations (|134p : 

Co = 2 {^uvuv ^u^)+-) — 2 i^-i K- Cuv + -} — Cy_^u — g {Ruvuv ~t" R-] h- ^~ ^ Rv-+u Ruv + -) , (140a) 

C_^_i — Cvuv- — C^-v- — hiRvuv- -\- R+-v-^ , (140b) 

UVU+ —+U+ UVU+ —+U+ 

C^a = Cv-v- = Rv-v- . (140c) 



U-\-U-\- + 
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Equations (|134p and (|140p together account for all the non-zero components of the Weyl tensor Ckimn in terms of Cs 
and its angular conjugate C*. 

The angular conjugates C* are proportional to spin harmonics of spin-weight —s, as follows from the fact that 
angular conjugation swaps angular indices + ^ — on any Riemann component, and the rule (j96p that the spin of any 
Riemann component equals the sum of the +'s and — 's of its covariant indices: 

c* (X . (141) 

Linearized expressions for the spin-s components Cs of the complexified Weyl tensor in a Newman-Pcnrose tetrad 
follow from their general expressions (|134p in terms of the Riemann tensor, and from the linearized expressions pi3p - 
(I117P for the Riemann tensor in terms of the tetrad connections and their directed derivatives. Because of the variety 
of expressions for the Riemann tensor, there are several ways to write the spin-0 and spin-±l components Co and 
C±i, though only one way to write the spin-±2 components C±2- 

Being tetrad-frame quantities, all components Cs of the complexified Weyl tensor (and their angular conjugates) 
are coordinate gauge- invariant, although the split between the unperturbed and perturbed parts of the polar part 
of Co, the only component whose unperturbed value is non- vanishing, is not coordinate gauge-invariant. The spin-0 
and spin-±2 components Co and C±2 (and their angular conjugates) are not only coordinate gauge-invariant but 
tetrad gauge-invariant with respect to all 6 arbitrary infinitesimal tetrad (Lorentz) transformations. The spin-±l 
components C±i (and their angular conjugates) are not tetrad gauge-invariant with respect to arbitrary infinitesimal 
tetrad transformations, being changed by a tetrad transformation that varies the antisymmetric part of ipv+ (or 

its angular conjugate (fiv-)- The following combinations of C±i or C*i and Co are tetrad gauge- invariant, but not 

u+ 

coordinate gauge- invariant: 

C±i + ^{(pv+ - ^+v)Co , (142a) 
Cli + Uvv^-V-v)Co ■ (142b) 

Special interest attaches to quantities that are both coordinate and tetrad gauge- invariant, and also vanishing 
in the background, because such quantities constitute physical perturbations that are unchanged by any gauge 
transformation. There are five components of the complexified Weyl tensor that are coordinate and tetrad gauge- 
invariant and vanishing in the background, namely the axial spin-0 component Cg"'' , 

Cq = ^(C'o ^ Cq ) = ^Cuv+- = \Ruv+- ) (143) 

and the spin-±2 components C±2 and their angular conjugates C*2. 

It follows from equations (|129p and (I142p that the following combinations of spin-±l Weyl and Einstein components 
are both coordinate and tetrad gauge- invariant, and vanishing in the background: 

K^uti + ^+-)^ ~ GyyGuu\C±i — SCoiGyu + C^ — )Gv+ + 3CoCtJtiCu+ , (144a) 
\{Gvu + C+_)^ — Gi,jjGuu'\C*_i — 3Cg(Gt,u + G^_)Cv- + 3CQCt>t)Gti- . (144b) 



B. Evolution of the Weyl tensor 



Equations governing the evolution of gravitational waves are obtained by applying the Bianchi identities to the Weyl 
tensor. The central result of this subsection is the set of linearized Weyl evolution equations (|156p . which describe the 
coupled evolution of the various spin s = — 2 to +2 components Cs of the complexified Weyl tensor in a spherically 
symmetric self-similar unperturbed background. These Weyl evolution equations can be combined in pairs to yield 
sourced wave equations (|165p for each of the spin components Cs ■ Since the different spin components are related to 
each other by the Weyl evolution equations p56p . the wave equations do not describe independent modes of vibration, 
but nevertheless it is useful to have the full suite of wave equations available. The next subsection, ^VC] considers 
the case of the spin-±2 components C±2, which describe the propagating components of gravitational waves, while 
the final subsections, W Gl and W HI consider the axial and polar spin-0 components Cq"'' and Cq^\ 

The Bianchi identities 



(145) 
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yield evolution equations for the Weyl tensor 

D^Cumn — Jlmn (146) 

with Weyl current 

Jlmn — \{DmGln — DnGiyn) — \{'^lnDmG ~ ^ImDnG) , (147) 

which satisfies the conservation law Ji„in — 0, in view of equation (|146p and the antisymmetry of Ckimn with 
respect to the indices kl. The corresponding evolution equations for the complexified Weyl tensor are 

D'^Cklmn = Jlmn (148) 

where Jimn is the complexified Weyl current 

Jlmn = ^ (S'^^S'^ + ^e™/"^ Jlrs , (149) 

which satisfies the conservation law D^Ji„in = 0. The complexified currents J/„m are related to the currents Jimn by 

Jvvu — Jv-+ — J-V+ — Tyijvvu Jv-+ \ — J~v+ J (150a) 

UUV u+- +M- ^ ^ IIUV tiH ' +U~ 

J+vu — J -\ h ~ Juv+ — ^ (j+vu ^" J^ h) ^ Juv+ 5 (150b) 

— uv — + - vu- -uv - + - vu — 

Jvv+ — Jvv+ 1 (150c) 
uu~ uu- 

J+V+ = J+V+ ■ (150d) 

—u— —u— 

Complexified Weyl currents Jimn with right bivector index mn equal to v— or u+ are zero 

J—v— — Juv— — Jvv— — t/^u— =0. (151) 
+U+ VU+ UU+ -u+ 

The angular conjugates J;*„„ of the complexified currents are given by the angular conjugates of equations ()150|) . 
which swaps angular indices + ^ — on the currents Jimn on the right hand sides (but does not swap angular indices 
on the complexified currents Jimn)- 

Jvvu — Jy-A- ~ J -v+ ~ ^(Jvvu (/d — + ) = J+v- , (152a) 

UUV u+- +u- «H — -u+ 

J+vu — J+ — \- — Juv+ — J-uv + J- + -) — Juv- , (152b) 

-uv — 1-- vu- +vu +-+ vu+ 

J VV+ ~ Jvv- , (152c) 

uu- UU+ 

JXv+ = J-v- . (152d) 

-u- +U+ 

In a general Newman-Penrose tetrad, the equations (|148p governing the evolution of the complexified Weyl tensor 
Cs are explicitly 

(153a) 
(153b) 
(153c) 
(153d) 

— a — <J —II — 

where the differential operators are 



2 r -\-vv 

-uu 






+ 0^v 

U 


ACo- 


iiA^^Cii - 


f r_„_ A^ C±2 


= A^ 


J -V + 

+u- 


2 r-i-ii+ 

-u- 


A^ 




+ oA± 


ACo- 


+1 A-u, A C+i - 

V 


+VV 


= A^ 


J + 

vu- 








-uu 


A^Co^ 


- ±1 Ai; A C±l 

u 


- ±2A^AC±2 


= A^ 


J VV-lf 

uu- 








3 r+i;+ 

— U- 


A^Co^ 


-±iA± ACii 


— +2^u A C+2 
v 


= A^ 


J+V + 

—u- 



Av = A^ 



±s 



dv + s Tuvv - s T-^v + (s + 3) r+„_ A-^ , (154a) 

±,A± = A' d^^s T-++ + s Tuv+ + (s + 3) r+^,„ A-^ , (154b) 
+- 
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which reduce in the unperturbed background to 



A — T> 



(155) 



Explicit expressions for the Weyl currents Jimn on the right hand sides of equations (jl53p are given in the Appendix, 
equations (jA2p . 

To hnear order of perturbations on the self-similar background, the Weyl evolution equations (|153p reduce (in a 
general gauge) to 



A- 



(dv + 3 r+, 

\ -uv / 







= A^ 


J-V + 

+u- 


— ±i'Du 

V 




= A^ 


J UV + 

vu- 




AC ^2 


= A' 


J vv-\- 
uu- 


— ±2T^u 


AC ^2 


= A^ 


J -iv + 



(156a) 

(156b) 
(156c) 
(156d) 



where the linearized Weyl currents Jimn on the right hand sides of equations (|156p are (again in a general gauge), 
from equations (|A2p . 



Jvvu — Jv- + 
uuv u+- 



J-V + 



= J+ 



f+vu — h 

—uv 



J uv-\ 
vu~ 

J vv^ 
uu- 



1 



dv {\GvU — (du ^ 2 T^VUU + Yj^u- — T-U+l GvV ~ A ^ q:l1^± Gv- ~\~ ^ ^ ±1^T Gv^ 

'II. ^ ' ^ V H.HV —Tij. Jj.y—' UU ii-L- 11 — 



+ (r-,,+ -r+._)(G„„ + G^_) 



dv (\Guu + \Gi\ A ^ ^:i'D± Gv- + r-u+ (Cyu + C I _) + V -u+ Gw 



d±{^Gvu + \ 



G, 



A ^^iDvGu+ — A ^^iDuGv+-\-A ^^2T^^:G++ 



(157a) 
(157b) 

(157c) 



-|- f r+tiii r + uv ) \Gyii H- j -|- r+u-u Gw Y' -\-vv Gv 

— uv —VU —VV UU —UU 1- 

1 1 

— ^± i'^Gyu H~ q ) "t" ^ tI^i' Gu+ — r+nu ( Gyu H" I — ) — T -^vv Guu -f- r_u+ C-u^- , (157d) 

' ^ U v~ -VU^ ' ~UU VV +v- u-j 

1 -| 

-i-i'Dv A. Gv+ {d± ~^ 2 'Vuv+ H~ r+uD ) Gvv — Y -i-vv [ Guu ~^ G^_ ) , (157e) 

U U- VU- -VU^ UU -UU ^ ^ J 

oVv A-^G++ - A-\iV^ Gv+ + r +f + (Gyu + G^_ ) + r+w-h Gvv . (157f) 



The equality of the expressions (|157ap and (jl57bp . and of (|157cp and (jl57dp . expresses the vanishing of the covariant 
derivative of the Einstein tensor, D"^Gmn — 0, which implies the conservation of energy-momentum. 
In the unperturbed background, the Weyl evolution equations (jl56p simplify to 



oVv A-^C = J.V+ , 

u +u- 

oV^A^^C = , 

where the only non-vanishing component of the unperturbed Weyl current is 





J vvu — J V — + — ~" 

uuv u+- 4 

1 

2 



Af)DvA ^ (^Gt,u + ^G+_) — A ^±2T^uGvv 



V UU 



J -v + 



(158a) 
(158b) 



\A-^'^^Vv A~^'^Gvu + \A-^qVv a G+_ + A-^Pu Gvv , (159) 

> U U V uuj ^ ' 

in which the equality of the two expressions expresses D"^Gmn = in the unperturbed background, which implies 



energy-momentum conservation in the unperturbed background. A "nice" expression for J- u+, obtained by combining 

+u- 

the first and second expressions of equations p59p in the proportions | and is 

1 r 



J -v+ — — 

6 



oI?ti A ^ (Cvu + Gj ) — ±2'Du A ^Gvv 

II ^ ' V 1IH 



(160) 
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In a general Newman-Penrose tetrad, the Weyl current conservation equations Z?' J;„m = are explicitly 

D Jluv = +l^u ^ Jvvu -l^v ^ Juuv A J-^vu + -lA^ A J-uv 

^ -uu Jvv+ ~t~ r_i_^^ Juu- ~t~ r_^_ tj^y^ r^^^ J -u- o ^ (i6ia) 

A? 1 , = 0^+ Jvvu — A? J+vu — ±2^M A? Jvv+ + ±2^L A? J+V+ — 2 Y +VV Juuv + 2 r+D+ J~uv = , (161b) 

where the differential operators o^m 



±s^v = A 



dv + s Tuvv - (s + 1) r-+„ + r-«+ + (s + 3) r+„- A-' , (162a) 

u VUU H u +u- -U+J 



a± - s r_++ + (s + 1) Tuv+ + T+uv + {s + 3) T+vu A"^ , (162b) 

H vu- -vu -uv] 

which reduce in the unperturbed background to 

,A:„ = ,1?™ . (163) 

Linearized, the Weyl current conservation equations (|16ip reduce (in a general gauge) to 

+iA'^ A^ Jvvu - -iKA^ Juuv- ^iT^-A^J^vu + -iV+A^ = 0, (164a) 
qA I A Jvvu A J+vu ±2^u A J vv+ ~\~ ±2^T A J +v+ 2 F+Dt; Juuv — . (164b) 

^ uuv u ^uv V uu- -ti- -uu vvu 

Combining the linearized Weyl evolution equations (|156p in pairs, by taking one of the operators of 
equations (|162p times one equation, minus another of the operators sA^ times an adjacent equation, the choice 
of operators being guided by the Weyl current conservation equations (|164p . yields equations that look like sourced 
wave equations for each of the (perturbations of the) spin components Cs of the complexified Weyl tensor (in a general 
gauge): 

\ u J v\ u +U-J ^\ VU-J 

f iipt" ^^Vu + iiP^^ ±iV^\ A (7±i = oAl f- oA. A-^C + A^ J_„+) - oVv {- oA± ^-^C + J„„+) , (165b) 

fiiX't"^^!?^ +^^pt± 2? "j^C* = ±22?«f3r+™C + A2 J^^+) -±22?:ff3r+„+C + A2 J+„+) , (165c) 

f±22?^"±22?n + ±2l?^=f ±2I?t')AC'±2 = ±il?± f 3 T C + J^^+ ) - ^3 T +„+ C + J+„+ ) . (165d) 

On the left hand side of equations (|165p . the leftmost of each pair of differential operators has been replaced with an 
equivalent Hcrmitian conjugate operator, per equations (j70p and (|62p . which makes manifest the Hermitian character 
of the wave operators on the left hand side. It should be emphasized that these are not independent wave equations: 
the various spin components Cs are not independently adjustable, but rather fluctuate in harmony with each other 
in accordance with the Weyl evolution equations (|156p . 

Equation (|165ap apparently constitutes two separate equations for Cq, while equations (|165bp and (|165cp apparently 
provide two separate equations for each of C±i. It is a straightforward if somewhat lengthy exercise to confirm 
directly, using the commutation relations (|32p. the operator relations (|55|) and ([75)) . and the Weyl current conservation 
equations (|164p . that the separate equations for each Cg are in fact equivalent. For Co, the confirmation also invokes 
the equivalence of the two expressions (lAlcP for while for C^i, it is necessary to use expression (jl34bp for C^i 

with the top lines of equations ( |AlgP and (jAlip substituted. 

It should be commented that it is legitimate to replace the differential operators ±iA^ and oA'_,_ on the right hand 

V 

sides of equations (jl65ap and l|165bp with their unperturbed limits ±iDu and qD^^ since the combined quantities in 
parentheses that they operate on are already of first order (although the separate terms inside the parentheses are 
not individually of first order). However, it is convenient to retain the operators in the unperturbed form, since this 
makes application of the Weyl current conservation equations (I164p . as in the previous paragraph, transparent. 

One of the more complicated steps in confirming that equations (I165bp and (|165cp are equivalent is to derive 
equation (|A5p given in the Appendix, which linearizes (in a general gauge) to 

foA'±oA„-oAUA±)A-iC = 3cf2ylC'±i-A±22?«A"^r+™ +±2V^T^v+\ - AT +vv o'Du A'^ C , (166) 

\ u u J V -uu -u-J -uu V 
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the last term on the right hand side of which can be converted to unperturbed Weyl currents Jimn using 
equation (|158ap . 

Another relation that proves useful in reducing the wave equation for the axial spin-0 component Cq^^ W GI below, 
is equation (jA6p in the Appendix, which linearizes (in a general gauge) to 



(+iA', oA, - +iA'_ oA+ - +iA:: oA: + +iA'* oA;) A-^C 



6CAC, 



(a) 



A-^C 



(167) 



the last set of terms on the right hand side of which can be converted to unperturbed Weyl currents Jimn using 
equation (|158ap . 



C. Spin-two (propagating) components of gravitational waves 

The spin-±2 components C±2 of the complexified Weyl tensor, and their angular conjugates (7*2, are of particular 
interest because they describe the ingoing and outgoing propagating components of gravitational waves. They 
are gauge-invariant with respect to infinitesimal coordinate and tetrad transformations, and vanishing in the 
background, and therefore represent physical perturbations. The main result of this subsection is the sourced spin- 
±2 wave equation (jl72p [and its angular conjugate (|178p ]. whose homogeneous solutions are separable as given by 
equations (|179p . the radial part satisfying the generalized spin-±2 Tcukolsky equation (|180p . 

The wave equation (|165dp governing the evolution of each of the spin-±2 components C±2 is, like C±2 itself, 
coordinate and tetrad gauge- invariant. However, there is an inherent ambiguity in deciding how to split the wave 
equation into a left and a right hand side - a "homogeneous" part and a "source" part - an ambiguity that is not 
resolved simply by writing the equation in a form such as equation (jl65dp . The problem is that, although C±2 
is itself coordinate and tetrad gauge-invariant, the space of perturbations "orthogonal" to C±2 is not. Indeed, as 
noted in ijlVEl the only coordinate and tetrad gauge-invariant perturbations to the Einstein tensor are its spin-±2 
components G++. For example, the Einstein component Gv+ [which appears in the Weyl currents Jvv+ and J+v+, 

uu- -u~ 

equations p57ep and (jl57f|) . on the right hand side of the wave equation (|165dp ] is not tetrad gauge-invariant, so by 
a suitable infinitesimal tetrad gauge transformation (that is, by a suitable choice of tetrad gauge) it is possible in 
effect to add to arbitrary amounts of +22^-^+2 (the angular lowering operator -^2^^- being required to convert the 
spin-2 object (5+2 into a spin-1 object like Gv+)- The word "arbitrary" in the previous sentence is a bit too general: 
it is desirable that the homogeneous part of the wave equation should be (coordinate and tetrad) gauge-invariant, 
so the amount of +22^-^+2 added to G„+ should be restricted to being gauge-invariant; but this still allows a broad 
range of choices. It would be possible to split Gy+ unambigously into homogeneous and source parts respectively 
proportional and orthogonal to C±2 if there were an unambiguous notion of orthogonality for C±2, which would 
be true if C-1-2 were an eigenmode of some operator. However, C±2, whose evolution is governed by the coupled 
Weyl evolution equations (jl56p . is not an eigenmode of those equations, except in the case of vanishing background 
energy-momentum tensor, the Schwarzschild geometry. 

The argument of the previous paragraph shows that splitting the wave equation for the spin-±2 components G±2 into 
homogeneous and source parts effectively involves choosing a gauge. Intriguingiy, spherical gauge, equations (|106p . 
which emerged in ? ijIVBI and IIV CI from the (not necessary, but nevertheless "natural" ) requirement that all tetrad 
connections be expandable in spin harmonics, provides a split that proves particularly elegant. In spherical gauge, 
the connections T+v+ and r-i-ii+ vanish, and C±2 and C*2 can be expressed in terms of tetrad connections as, 

-u- -~v~ 

equations P^ic)) and (fHOc)) with (|117bp and pTfc)) . 

C±2 = A-\iV±r+uu , (168a) 

-UU 

CI2 - A-^^{D^Y-uu . (168b) 

For subsequent reference it is useful to record here that in spherical gauge the spin-±2 components G++ of the Einstein 
tensor are, equations (jl25f|l and (|117ap . 

G++ = -2A-i±iI?±r+™ . (169) 

— -uv 

Although the right hand sides of equations (|168ap . (jl68bp . and (|169p are not (tetrad) gauge-invariant, they 
are evidently equivalent to gauge- invariant quantities. For non-dipole spherical harmonics, with I ^ 1, 
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equations (|168ap . (jl68bp . and (|169p invert to yield expressions for T+vv, T^w, and r+i,„ (in spherical gauge) in 

-UU +UU -uv 

terms of C±2, C22i 

2A 
2A 

The right hand side of the wave equation (jl65dp for C±2 contains the connections T +vv in two places, in a form that 

-uu 

converts them precisely into C±2, and it makes sense to take these terms over to the left hand side of the wave equation, 
where they become incorporated into the wave operator. The first of the two terms is proportional to the Weyl scalar 



C, while the second occurs in the current Jw+j equation (1157ep . and is proportional to i5 = ^A^iG^u + The 
two terms combine in the proportions 

- iC + E . (171) 
The resulting wave equation for the spin-±2 Weyl tensor C±2 is 

(±22?^" ±2Vu + ±22?^^ ±2V^ -iC + E)A C±2 = 5±2 , (172) 

where the spin-±2 source term 5*^2 on the right hand side is, in spherical gauge, 
A r 



A-^ ±iVv A^ oVv A-^G++ + 2 ±iVv ±iX>± G„+ - A ±iX>± (a± + 2 r„„+ + r+™)G'™ . (173) 



Although spherical gauge has been invoked to split the wave equation (jl72p into homogeneous and source parts, the 
left hand (homogeneous) side of equation (jl72p is coordinate and tetrad gauge-invariant, and therefore the spin-±2 
source term £"±2 constituting the right hand side must be equivalent to a coordinate and tetrad gauge-invariant object. 
In fact the full (gauge-invariant) expression for the source term S±2 is (in a general gauge, not spherical gauge) 



'±2 — — ±lT^v A'^ \T (Gyu -\- G I _) -|- i r+ji+ Gt;t)] -|- i A'^ r+ii+ ±2^" ^ ^Gvv 

u L -u- ^ I / z -V- UU-' z -u- V uu 

A 
+ 2 



- A-^ ±{Dv A^ oVv A-^G++ +2±iVv ±iV±Gv+ - A±iVJd± +2Tuv+ +T+uv)Gvv . (174) 

1A II. 7/, II _ ^ _T1U' till. 



It is evident from equation (|174p that spherical gauge, where r+^+ and T+u+ vanish, indeed leads to a significant 

-U- -V- 

simplification of the source term. 

Spherical gauge, equation (|106p . does not fix the gauge uniquely, since it still leaves two infinitesimal gauge freedoms, 
corresponding to 1 coordinate freedom in the choice of radial coordinate r, and 1 tetrad gauge freedom in the 
antisymmetric part of ifyu (a Lorentz boost in the radial direction). It is worth pointing out that these gauge 
freedoms are precisely sufficient to set the contribution from the unperturbed part of Gvv in the source term 5'±2, 
equation (|173p . to zero (although in the next paragraph it will be argued that this may not be a preferred gauge 
choice). That is, the contribution to the source term S+2 from Gw splits into parts depending on the unperturbed 

uu 

1 

Einstein component Gw and its perturbation G 

UU UU 

du + 2 H~ -^uvjGvv ~h -A 0^-|- Gvv . (175) 



Since one of ip+u or ip+v can be adjusted with the coordinate gauge freedom in e" or e", equations (|102ep and (|102f|l . 
while one of Tuv+ can be adjusted with the tetrad gauge freedom in Lp^u, equation (|104hp . the two gauge freedoms 

vu- 

available in spherical gauge are just such as to make it possible to set 



[ij^+udv + (p+vdu + 2r™+ -I- r+uv)Gvv — . (176) 

^ -V u ' -u V vu- ~vu' uu 
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If this gauge choice is made, then the spin-±2 source term becomes 



5* 



±2 



A-^±iVv oVv A-^G++ + 2 ^{Dv ±iV± Gv+ - ±iX>± qV+ 



(177) 



which is precisely the form that the spin-±2 source term takes in the case of the Schwarzschild geometry, where ah 
unperturbed components of the Einstein tensor vanish, and all perturbations to the tetrad frame Einstein tensor are 
coordinate and tetrad gauge-invariant. 

However, the gauge conditions (|176|) . though feasible, are ad hoc, unlike the spherical gauge conditions (|106p . 
which have a "natural" origin (that all tetrad connections be expandable in spin harmonics). Moreover the gauge 
conditions (jl76p represent a rather complicated condition on the vierbein perturbations v?mrij which may not be ideal 
for numerical calculations. Indeed, in a hypothetical asymptotically flat empty region far from the black hole, the 



unperturbed Einstein component Gw vanishes asymptotically, so that the conditions (|176p are satisfied automatically 
rather than being gauge conditions, and then how the conditions (jl76p translate into conditions on the vierbein 
perturbations Lpmn may become a numerically delicate issue. 

Equation (|173p shows that the tensor (spin-±2) components of gravitational waves couple not only to tensor sources 
(G++) but also to vector {Gv+) and scalar (Gw) sources. If there is a non-vanishing background energy-momentum, 

— u— uu 

then gravitational waves can excite other (non-gravitational) kinds of waves in the energy-momentum, which in turn 
provide a source for gravitational waves. The coupling of modes in black holes is unlike the case of the Friedmann- 
Robertson- Walker metric of cosmology, where spatial translation and rotation symmetries ensure that scalar, vector, 
and tensor modes decouple completely from each other [l|, 0, H^l • 

The angular conjugates (7*2 of the spin-±2 components of the complexified Weyl tensor satisfy the angular conjugate 
of equation (|172p 



±21 



Vu + ^2V^^ ^2V± -iC + E\A CI, = SI 



=f2 



±2 



^±2 



(178) 



The angular Hermitian operators +22?^" +22^- and -22^^^-225+ inside the wave operators of equation (|172p and its 
angular conjugate (|178p have the same eigenvalues (with angular eigenfunctions that are complex conjugates of each 
other), in view of equations (j63p . and it follows that the temporal-radial (i.e. non-angular) part of the eigenfunctions 
of the wave operators are the same for both C±2 and its angular conjugate C*2- However, the temporal-radial parts 
of the source terms S^2 are in general not equal to their angular conjugates >S'*2- 

The wave operator on the left hand side of equation (|172p differs from the wave operator discussed in ijlll El onlv by 
the addition of the terms — ZC + E, which are functions only of conformal radius r. The homogeneous solutions (when 
the source terms vanish) of the wave equation (|172p and its angular conjugate (|178p are, similarly to equation (|78p . 



C±2 ^ = (^7e")TiA-3e— 



±2 



±2^1 



(179) 



in which C±2 and its angular conjugate differ only in that their angular parts are complex conjugates of each other. 
Note that the complex conjugate of a harmonic of spin s has opposite spin — s, in accordance with equation (j68p and in 
agreement with equation (|14ip . The polar and axial parts of the eigenfunctions (|179p have the same temporal-radial 
dependence, and simply project out the real and imaginary parts of the angular factor, the spin-weighted spherical 



harmonic ±2^^;?! 



The radial function ±2Cu}i{i') in equation (jl79p depends only on conformal radius r, and satisfies 



the second order ordinary differential equation 

92 



dr 



,*2 



iuj ± 2n) + [/(/ + l)-4-18C + 2£;]A^±2C^/ =0 , 



(180) 



which agrees with the generalized Teukolsky equation ^ given in the Introduction, since, equations (jH]) and p27p . 

- 18C -I- 2£; = -5(i? - P) - 4Pl + 18M . (181) 



Equation (jl80p has two linearly independent eigensolutions, which may be identified as ingoing and outgoing based 
on their properties in a hypothetical asymptotically flat region at large radius, as already discussed in ijlllFl In a 
hypothetical aysmptotically flat region at large radius, the radial eigensolutions sCu)i (with s = ±2) go over to the 
asymptotic solutions gZ; discussed in ijlllFl equation 



{l + s)l 



1/2 



1/2 



{2Lor*ye- 



as r 



00 



(182) 
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Because the wave operators for modes of opposite spin are complex conjugates of each other, equation (|180p . it follows 
that the ingoing and outgoing radial eigensolutions sC^i of opposite spin can be taken to be complex conjugates of 
each other, in complete analogy to equation ([86|l. 

s~Ci^{-scZTT ■ (183) 

The propagating components of Cs (x e~^'^* sCuii are those that fall off most slowly at large radius, which identifies 
the positive spin component C+2 as describing the propagating component of ingoing gravitational waves, and the 
negative spin component C_2 as describing the propagating component of outgoing gravitational waves: 

C% cx A-^e-~''^\2iZ - A-3(r*)2e-*'^(*+'-') (propagating, ingoing) , (184a) 
(7°^* (X A-3e-*'^*_2C* - A"3(^*)2g^».^(t-r-) (propagating, outgoing) . (184b) 

The other of the two modes for each spin corresponds to the non-propagating, short-range spin-(— s) partner of a 
propagating mode of spin s: 

C% « A-3e-'"*_2C; - A-3(r.*)-2e-»'^(*+'-*) (non-propagating, ingoing) , (185a) 
Clf oc A-^e-"^\2C°^' A-3(r*)-2e-*'^(*-r-*) (non-propagating, outgoing) . (185b) 



D. From spin-±2 to general spin-s components of gravitational waves 

The Weyl evolution equations (|156p relate the five spin s = —2 to +2 components Cs of the complexified Weyl tensor 
to each other, so that once one component has been specified, then the others follow. Thus, except for monopole and 
dipole harmonics, Z = or 1, for which all spin-±2 quantities vanish identically, and which will be considered specially 
in jji jVEl and IV Fl either of the spin-±2 components C^2 of the complexified Weyl tensor is sufficient to determine 
all the other Weyl components. Starting with either of the components C+2 or (7-2, the Weyl evolution equations, 
reordered in the sequence (|156d|) . (|156b|) . (|156ap . and (|156cp . yield in succession 

C±2 ^ C±i (5o ^ C^i ^ C^2 ■ (186) 

The angular raising/lowering operators sT^± involved in extracting each successive spin-s component can be canceled 

with corresponding angular lowering/raising operators s+i2?^, and then equation (jl56dp yields C±i in terms of C±2 

1 

(plus other things), equation (jl56bp yields Co in terms of C±i (plus other things), and so on. 

It is not worth writing out explicitly all of the resulting set of expressions for the various spin-s components Cs 
of the Weyl tensor, because the expressions for the most part offer little insight beyond that already in the Weyl 
evolution equations (|156p . 

However, the first rung on the ladder, equation (|156dp . which yields C±i in terms of C±2 plus other things, 
does deserve comment, because it reveals explicitly how the spin-±l and spin-±2 components of the Weyl tensor 
fluctuate together rather than varying independently. In spherical gauge, and for non-dipole modes, the Weyl evolution 
equation (jl56dp integrates to 

C±i + (I- r){l + 2) ^'^^ (-^''±2Pn A C^2 + \AoV.^ [l^l), (187) 

which expresses the combination C^i -f hGv+ of spin-±l components of the Weyl and Einstein tensors in terms of the 
spin-±2 components C±2 and G++ of the Weyl and Einstein tensors. Equation (|187p thus shows that, as asserted, the 

spin-±l and spin-±2 components of the Weyl tensor fluctuate together rather than varying independently. In spherical 
gauge, the spin-±2 components of the Weyl and Einstein tensors are related to tetrad connections by equations (|168p 
and (I169p . and equation ()187p can be rewritten as 

C±l + \Gv+ — A ^\ ±2'Du T -f-vv — qT^v T+vu I , (188) 

which can be recognized as the same as the Riemann equation (jll6ap . expressed in spherical gauge (where V^m = 
r+uu), and is therefore valid for dipole modes / = 1 as well as non-dipole modes. 

— vu 
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E. Gravitational monopole modes 

For monopole or dipole modes, ^ = or 1, all spin-±2 quantities vanish identically, and therefore the spin-±2 
components C^2 of the c omp lexified Weyl tensor cannot describe these modes. Thus monopole and dipole modes, 
which as emphasized by [9l| are essential to a complete description of gravitational perturbations, require special 
attention. The present subsection covers monopole modes, while the following subsection, ijVFl will deal with dipole 
modes. The main results of this subsection are: the polar monopole equation (|193p . which simplifies in spherical 
gauge to equation (|194|) ; and the axial monopole equation ()195|) , which states simply that the axial monopole vanishes 
identically. 

An equation for the monopole follows from the first Weyl evolution equation (|f 56a[) . After some manipulation, the 
Weyl evolution equation (|156a|) can be recast in the following form (in a general gauge, and for arbitrary harmonics), 
with the feature that every term is individually of first order: 

+ r_„+ + r+„_)A-2AA + oi?«yic'^''^ - ^{D^AC^x = i(r_„+ - r+„_)(3C + £;) - iA2(r+„_ -r_„+)G™ 
+ i/3«(±ii?^r+™ +^ii?±r_™) + i/3,;(±ii?^r+™ + ^ii?±r_™) - \a^{d^g-„- , (i89) 

where Cq""* is the axial spin-0 component of the Weyl tensor, and M is defined to be the following dimensionless 
combination of polar spin-0 components of the Weyl tensor, Einstein tensor, and tetrad connections 

M = A^ {c^^^ - ia„ - iG+_ + \ (r_,^ + r+„_) (r+„_ + r_„+)] . (i90) 

The combination of polar spin-0 components of the Weyl and Einstein tensors in the definition p90p of M is 

proportional to one of the Riemann components, A^{c'^^f^ — \Gvu — — \J^R ^-^- , which in the unperturbed 

background reduces to minus the dimensionless interior mass, — Af, equation (jll8b[) . The polar spin-0 quantity M 
is just iv4^ times the object on the left hand side of equation (|f 24p . so the expression for M in terms of vierbein 
perturbations (/?„in is given by times the right hand side of equation (jl24[) . The quantity M is not coordinate 
gauge- invariant, although A~'^M is coordinate gauge-invariant, and A/" is not tetrad gauge-invariant, being changed by 
infinitesimal tetrad gauge transformations of any of the antisymmetric parts of 93+1,, V'+u, H^-v, or (/3_u, equation (|f 24p . 
However, the monopole part of A/" is tetrad gauge-invariant, since iy9+u, <~P-v^ and 93_„ are spin-±f objects, whose 
monopole parts therefore vanish identically. The notable asset of M is that the split between its unperturbed and 
perturbed parts is coordinate and tetrad gauge-invariant, the unperturbed part being simply a constant, a coordinate 
and tetrad gauge-invariant object, 

k = -\. (i9f) 

For the monopole mode, I = 0, all non-zero spin objects in equation (|f 89p vanish identically [including V — V^y_ 
and r^„_ — r_u+, whose expressions (If 04cp in terms of vierbein components <^ran contain only spin-±f terms], and 
the equation turns into one that determines the monopole to linear order: 

A^fa. -)-r_.+ -)-r+„_)A-2_^ + oP«AC'^"' = (; = o), (f92) 

a coordinate and tetrad gauge-invariant equation. 

The polar part of the monopole equation p92p is, in a general gauge, 

^2 fa„ + r-„+ + A-^AA^ (; = o). (f93) 

In spherical gauge, where equation (|f f 2ap holds, equation (|f 93p is equivalent to dvM = 0, which integrates to 

AA=-i (; = 0), (f94) 

the constant of integration being determined by the unperturbed value (|f 9f p . The fact that in spherical gauge the 
monopole component of the polar spin-0 quantity Af equals — ^ not only to unperturbed order but also to linear order 
is a pretty result, and definitely another feather in the cap of spherical gauge. This feature will be exploited in W Ifl 
below to derive a wave equation for polar spin-0 modes. 

The axial part of the monopole equation (|f92p is, in a general gauge, qT^v ACq"'^ — 0, which integrates immediately 

to 

C^'^) =0 {1 = 0) , (f95) 

that is, the monopole axial spin-0 component of the Weyl tensor is identically zero, a coordinate and tetrad gauge- 
invariant statement. 
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F. Gravitational dipole modes 

An equation for the dipole follows from the third Weyl evolution equation p56cp with (7-1-2 set to zero, which gives 
(in a general gauge) 

iiVv A Cii = (3C - E) r + ™ + 1^2 ^{Dv A-^Gv+ -\A^{d± + 2 Tuv+ + T+m^Gw {l^l) ■ (196) 

In contrast to equation ()188p . which constitutes an identity between Riemann components and tetrad connections, 
equation (jl96p is a genuine evolutionary equation, a first order differential equation governing the dipole spin-±l 
components C±i of the Weyl tensor. It is tempting to think that the tetrad connections T +vv in the dipole evolution 

~uu 

equation (|196p would vanish in spherical gauge, in view of the relations (|170ap between these tetrad connections and 
the spin-±2 components C±2 of the Weyl tensor, but for dipole modes the relations (|170a|) fails, and so the tetrad 
connections do not necessarily vanish. 

G. Axial spin-zero component of gravitational waves 

The previous four subsections, CtfVFl have given a complete description of gravitational waves (perturbations 
to the Weyl tensor) in spherically symmetric self-similar black holes, including monopole and dipole modes, so it 
might seem superfluous to explore the equations any further. 

However, there are two good reasons not to stop at this point, but rather to consider wave equations for the 
polar and axial parts of the spin-0 component Co of the complexified Weyl tensor. The first good reason is that the 
potentials in the wave operators for the polar and axial spin-0 components are real, not complex. For real potentials, 
theorems familiar from quantum mechanics apply, such as whether or not there exists a discrete spectrum of bound 
states (which happens if the potential has one or more minima). And the second reason is that it is desirable to 
attempt to make contact with the traditional description of perturbations in terms of polar and axial perturbations 
and the associated Regge- Wheeler and Zerilli potentials [H, IH, [H, [ol, . 

As matters stand, the connection between the Newman- Penrose formalism and the traditional polar-axial description 
of perturbations remains incompletely understood. Chandrasekhar [l8l |. by an ingenious suite of manipulations, is 
able to work from the spin-±2 wave equations in the Newman-Penrose formalism to derive a non-linear differential 
equation, his equation (311) on page 186, which is satisfied by the polar and axial potentials of the Schwarzschild and 
Reissner-Nordstrom geometries, and which is such as to ensure that the potentials yield the same transmission and 
reflection coefficients. Chandrasekhar's primary motivation was to recover the Regge- Wheeler and Zerilli decoupled 
wave equations for the polar and axial perturbations of the Schwarzschild and Reissner-Nordstrom geometries, and 
in this he was successful. However, it seems unlikely that the gravito-electromagnetic wave equations would decouple 
in the case of general spherically symmetric self-similar solutions, and it is unclear how Chandrasekhar's procedure 
might provide a way forward in this more general case. 

On the other hand, the Newman-Penrose formalism does volunteer a spin-0 wave equation, namely equation (|156ap . 
and it would seem good practice to explore where this equation leads. Chandrasekhar p^ . [19] did not consider this 
spin-0 wave equation. 

For non-zero spin s, the two distinct degrees of freedom of a Weyl component Ca are represented by ingoing and 
outgoing modes, which satisfy wave equations with complex conjugate potentials, and the potentials are the same for 
both polar and axial modes. For zero spin, by contrast, the two distinct degrees of freedom of Cq are represented 
by polar and axial modes, which satisfy wave equations with two different real potentials, and the potentials are the 
same for both ingoing and outgoing modes. The present subsection derives the wave equation p99p governing the 

axial spin-0 component Cg"-* of the Weyl tensor, a much easier task than for the polar spin-0 component c'^^ , which 
is considered in the next subsection, WHI 

The reason that the axial spin-0 component Cq"' , equation (|143p , is straightforward to address is that it is coordinate 
and tetrad gauge-invariant and vanishing in the background, the only other component of the Weyl tensor that has 
this happy property besides the spin-±2 components C±2 and C*2 considered in W CI A wave equation for the axial 

spin-0 component Cg"' follows from the wave equation p65ap minus its angular conjugate. 

Somewhat surprisingly, spherical gauge, which proved to be a reliable guide in the case of the spin-±2 wave equations, 
turns out not to provide the most elegant gauge for splitting the axial spin-0 wave equation into homogeneous and 
source parts. In fact, as will be seen below, there is a gauge in which the axial spin-0 source term takes its (apparently) 
simplest form, but not only is this gauge not spherical gauge, but it is incompatible with spherical gauge. The gauge, 
which can be termed "axial gauge" , consists of the two conditions 



(197) 
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which can be accomplished for example by 2 infinitesimal tetrad transformations of (f^y and (f-u, in accordance 
with equations (|104cp . As discussed in jjlVDi the requirement that perturbations be expandable in spin- weighted 
harmonics of itself imposes no gauge conditions, so the axial gauge conditions (|197p are entirely legitimate. But it is 
curious that the spin-±2 and axial spin-0 wave equations should "prefer" two mutually incompatible gauges. In the 
next subsection, ijVHl it will be found that the polar spin-0 wave equation "prefers" spherical gauge. 

In carrying out the somewhat lengthy calculation needed to evaluate the source term for the axial spin-0 wave 
equation, one finds that (much as occurred earlier for spin-±2 waves) there are two sets of terms that naturally 

convert themselves into quantities proportional to C'q"'' , and which are therefore incorporated naturally into the wave 
operator on the left hand side. One of these two sets of terms, proportional to the unperturbed Weyl scalar C, yields 

Cq via equation p67p . The other set of terms arises from the terms proportional to £' = ^A^l^Gyu + G+_) in 
the Weyl currents J-v+ and Juv+, equations (|157bl) and (I157dll . In axial gauge, the terms prove to combine in the 
proportions 

3C -E , (198) 

which is opposite in sign to the spin-±2 expression (|17ip . 

The resulting wave equation for the axial spin-0 component Cq"' of the complexified Weyl tensor is 

(o2?^" oVv + oV^^ o2?± +3C-E^A C^"^ = 5^"' , (199) 

where the axial spin-0 source term Sq"""^ is, in axial gauge, 

(a) ^ 

Equation (|200p is precisely the form that the axial spin-0 source term takes in the case of the Schwarzschild geometry, 
where all unperturbed components of the Einstein tensor vanish, and all perturbations to the tetrad-frame Einstein 
tensor are coordinate and tetrad gauge-invariant. The full (gauge-invariant) expression for the source term Sq""^ is (in 
a general gauge, not axial gauge) 

'S'q — — |+ii?u A [(r_i,^ — r^.„„) (Ciju + G_i__) — (r^„_ — r_„_i_)Gt,„] 



+ -iT^v ^[(r+u- — r„„^) (g^ju 4- G_|__) — (r_i,+ — ^+v-)Guu] 

+ +iDu Gy+ — +iDu ^{D^ Gy_ + ^iDy Gu- — -iDy +iD_ Gu+ 1 , (201) 

and it is evident that, as asserted above, the axial spin-0 source term 5*0"' takes its apparently simplest form in axial 
gauge, equations (|197p . 

The special case of the monopole mode, I = 0, was considered already in S|VEl where it was found that the axial 
monopole mode vanishes identically, equation (|195p . For the monopole mode, the axial gauge conditions (|197p are 
automatically satisfied, and do not constitute additional gauge conditions, since r_^_,_ — T^y__ and T^u- — ^-u+ depend 
only on spin-±l components of the vierbein perturbations ifmn, equations (|104cp . and therefore automatically vanish 
identically for the monopole. 

The wave operator on the left hand side of equation (|199p differs from the wave operator discussed in ijlllEl onlv by 
the addition of the terms 3G — E, which are functions only of conformal radius The homogeneous solutions (when 
the source terms vanish) are, similarly to equation (j78p . 

g'"^' = A-^e-^-\£l^^ (r) oY,m {0, <P) . (202) 

Previously, in the case of the spin-±2 components G±2 of the Weyl tensor, the polar and axial parts of the 
eigenf unctions (|179p projected out the real and imaginary parts of the spin-±2 spherical harmonic ±2Yim- The 
two parts could be distinguished fundamentally because the complex conjugate of a harmonic of spin s has opposite 
spin — s, equation (|68p. In the present case of the axial spin-0 component Cq°'\ the complex conjugate of a spin-0 
spherical harmonic is just another spin-0 spherical harmonic, so the real and imaginary parts are not fundamentally 
distinguishable. Thus the axial spin-0 eigenfunctions (I202p encompass both real and imaginary parts of the spin-0 
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spherical harmonic o^Zm, not just the imaginary part. The radial function oc^Si{r) in equation (|202p depends only on 
conformal radius r, and satisfies the second order ordinary differential equation 

-Lo^+ [l{l + 1) + 6C - 2£;] A I o^L'',^ = . (203) 



The effective potential + 1) + 6C — 2£'] A of the wave operator for qc^";'' is real, in contrast to the effective potential 
of the wave operator for ±2Ccoi, equation (|180p . which is complex. In terms of the dimensionless interior mass M, 
equation and dimensionless proper density R and radial pressure P, equations ([1^ . the factor 6C — 2E in the 
potential is, equations ((44l) and (|127p . 

6C -2E = R - P -6M . (204) 

Equation (|203p has two linearly independent eigensolutions, which may be taken to be complex conjugates of each 
other, and which may be identified as ingoing and outgoing based on their properties in a hypothetical asymptotically 
flat region at large radius, as discussed in mil Fl In a hypothetical aysmptotically flat region at large radius, the 
radial eigensolutions qc^"' go over to the asymptotic solutions qZi discussed in i jlllFl equation ([57)l . 

-.(a) 111 in( *\ p-iu/r* 

o'^uji ^ oZi [ujr ) ^ e ^* ^ (205) 

In the case of the Schwarzschild geometry, the potential in the wave operator of the wave equation (|203p for the 
axial spin-0 modes is precisely the familiar Regge- Wheeler [l^ [s^l potential V''^\ 

+ 1) + 6C - 2^;] A = [l{l + 1) - 6M] (1 - 2M) = V^*"' , (206) 

where the dimensionless interior mass M of the Schwarzschild black hole is related to its constant mass m, by 
M = m,/A ~ m,/r. The axial potential + 1) + 6C — 2i?]A also agrees with that obtained by [their 
equation (149)] for spherical stationary ideal-fluid imcharged "stars", and by [i^ [his equations (91)"(93) for s = 2] 
for general spherical stationary distributions. 

In the case of the Reissner-Nordstrom geometry, the axial potential + 1) + 6C — 2E\ A in the wave operator of 
equation (|203p is not the same as the Zerilli axial potential V^~^ p^.[63l[93t. This is because the wave equation (|203p 
does not describe decoupled eigenmodes of gravito-electromagnetic waves, but rather it describes homogeneous 
solutions of a wave equation (|199p whose source s'^^ includes electromagnetic perturbations that do not vanish 
for the Reissner-Nordstrom geometry. Physically, if the black hole is charged, then electromagnetic waves constitute a 
source for gravitational waves, and correspondingly gravitational waves constitute a source for electromagnetic waves 
[as will become evident below when electromagnetic waves are treated, for example equation (|270p ]. In order to 
obtain a decoupled wave equation, it would be necessary to flnd axial eigensolutions of the coupled wave equations for 
gravitational waves and electromagnetic waves. Thanks to the work of Zerilli and Moncrief 63., ,9^ it is known that 
such a decoupling can be accomplished for the Reissner-Nordstrom geometry, but for the case of general self-similar 
solutions there is no reason to expect that such a decoupling would occur. The present paper does not pursue this 
issue further. 



H. Polar spin-zero component of gravitational waves 



The derivation of a wave equation for the polar spin-0 component of gravitational waves presents a formidable 
challenge. This is largely because the polar spin-0 component Cq of the Weyl tensor does not vanish in the 
unperturbed background (even in the Schwarzschild geometry) , and consequently issues of gauge freedom are especially 
delicate. Nevertheless, for the reasons described at the beginning of the previous subsection, ^\ G[ it seems wise to 
pursue the polar part of the spin-0 wave equation (|165ap that arises in the Newman-Penrose formalism, to see where 
it leads. Where it leads is, apparently, to the sourced polar spin-0 wave equation (|217p . The present paper does not 
address the problem of how this equation and the other equations to which it couples might be decoupled so as to 
recover the Zerilli-Moncricf equations [H, in the Reissner-Nordstrom geometry. 

The polar spin-0 component Cq^-* = ■^(C'o + C'o) of the complexified Weyl tensor is gauge- invariant with respect 
to infinitesimal coordinate and tetrad transformations, but the split between the unperturbed and perturbed parts 
is not coordinate gauge- invariant. One might think (for an instant) that, since the unperturbed part of c'^^ is pure 
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monopole (1 = 0), at least the non-monopole parts of the expansion of Cq^-* in spherical harmonics would be gauge- 
invariant, but this is not true. Under an infinitesimal coordinate gauge transformation —^ x'^ = x^ + Ae^, the 



unperturbed part Cq of the unperturbed spin-0 Weyl tensor transforms as 

h^^\t, r) ^ hlj'\t', r') = dl^^ + Aie^d. + e"9„)4^) , (207) 

1 1 

so that, since Cq^'' = Cq^-* + Cq^^ as a whole is gauge-invariant, the perturbation C'q'''' transforms as 

^ ^(P) _ Aie"d, + e^du0P^ . (208) 

The dimensionless coordinate perturbations and e" in the transformation pOSp will in general contain arbitrary 

1 

non-monopole harmonics, changing the spherical harmonic expansion of the perturbation C'Jf'' arbitrarily. Indeed, 



as long as the unperturbed spin-0 Weyl tensor Cq is non-constant, which is generally true except in a hypothetical 

asymptotically flat region at large radius, it is possible in principle to eliminate the perturbation Cq to the polar 
spin-0 Weyl tensor altogether by a coordinate gauge transformation arranged to satisfy 

C^P^ - A(e''a„ + e"a„)C^^^ = . (209) 

However, it is not possible to eliminate simultaneously all polar spin-0 perturbations. Specifically, it is not possible 
to set the polar spin-0 vierbein perturbations fw, Vuu, and ip^u + Vuv simultaneously to zero, because each of 
these 3 vierbein perturbations varies only under a coordinate transformation of e" or e" (not under a coordinate 
transformation of e+ or e", nor under a tetrad transformation), equations p02p . and 3 gauge conditions cannot be 
achieved with just 2 freedoms. Thus there is reason to expect that the polar part of the spin-0 wave equation p65ap 
is not vacuous, but may in fact encode the evolution of polar spin-0 perturbations. 

A clue as to how to proceed is provided by the monopole mode of the polar spin-0 perturbation, which was found in 
SjVE]to be described most naturally not by the polar spin-0 Weyl component Cq''"' itself, but rather by the related polar 
spin-0 quantity M defined by equation (|190p . Although M is itself neither coordinate nor tetrad gauge-invariant, the 
equation (jl93[) governing the monopole part of A/" is coordinate and tetrad gauge-invariant. Therefore, if the monopole 
equation (I193P is subtracted from the general equation ()189p governing A/", then the result, being a difference of 
two coordinate and tetrad gauge-invariant equations, will be coordinate and tetrad gauge- invariant, and moreover 
vanishing in the background [that is to say, equation (jl89p as a whole is coordinate and tetrad gauge- invariant, though 
the split between its left and right hand sides is in general not coordinate and tetrad gauge- invariant] . 

As was found in §VEI the monopole part of M takes a particularly simple form in spherical gauge, being equal 
to its unperturbed constant value — ^ not only to unperturbed order but also to linear order, equation (jl94p . Thus 

in spherical gauge, subtracting the monopole equation from the general equation governing M amounts to replacing 

1 1 
A/" by its perturbation A/", a trivial operation. Although the resulting equation for the perturbation M is no longer 

gauge-invariant (being expressed in a particular gauge, spherical gauge), nevertheless the equation is equivalent to a 

1 

gauge-invariant equation, and therefore has real physical significance independent of any gauge. The perturbation M 
is coordinate and tetrad gauge-invariant with respect to the two gauge freedoms left in spherical gauge. 

The two Weyl evolution equations that combine to produce the spin-0 wave equation (|165ap are equations (|156ap 
and (jl56bp . The first of Weyl evolution equation (|156ap was already recast in terms of N in equation (|189p . Similarly 
recasting the second Weyl evolution equation (|156bp in terms of N gives (in a general gauge, and for arbitrary 
harmonics) 

9±AA + qI?± AC^"' -±iI?„AC±i ^ -A-i9±A-3r+.™C - iA2(a± -|-r+„„)(C„„ + G+_) 

+ iA2 5±[(r_.+ + r+ V-^ + U- -|-r-li+)l r^,^;!) GuU -\- h Zfl^V Gu+ -\- h A [5u Gv+ ■ (210) 

The monopole part of equation (|210p vanishes identically, since all the terms of the equation are spin-±l, so there is 
no need to subtract the monopole part of equation (|210p . 

As discussed above, it is helpful from this point to work in spherical gauge, equations (jl06p . because the subtraction 

of the monopole is trivial in this gauge, amounting to a replacement of J\f by its perturbation TV. In spherical gauge, 
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the perturbation A/" is, from ^A^ times equation (|124p . 

A/" = i [+lV^ 'P+v + Pv ^+u) + {Pu f-v + Pv f^u)] , (211) 



which may also be written 



Af ^ ^{+iV_d^A+ _iV^d_A) . (212) 



In spherical gauge, the term on the right hand side of equation (|210p involving the directed angular derivative of 
products of tetrad connections can be rewritten, using equations (|112ap and the commutation relations (j32[) . 



A2a±[(r-„+ +r+.-)(r+„- + r-„+)] = {(3^d, + Mu)d±A- /s^r^™ - /3,'r+„„ - 2/3„/3„r+™ . (213) 



1 /i2 
4 



When equations p89p and l|210p are combined into a polar wave equation for the perturbation A^, terms proportional 

to dj-A, namely the term —A~^dj-A on the right hand side of equation (|210p . and the terms proportional to d±A 

1 ^ 

on the right hand hand side of equation (|213p . convert into quantities proportional to Af, equation (|212p . and it is 

natural to incorporate these terms into the wave operator on the left hand side. 

1 

The resulting wave equation for the perturbation M is 

(+il?„ + - +iV_ + a_ + 2 Mv + 2 Pydu - 2A-i) = 25*^^^ , (214) 

where the factor of 2 on the right hand side is introduced so that it disappears from the wave equation (|217p below. 
The polar spin-0 source term S'g^^ on the right hand side of equation (j214p is, in spherical gauge, the unprepossessingly 
complicated object 

"So"^ = l{A^{+iV_d^ + _iV^d_){G,u + G^_) 

+ A-^ (/3„ o2?« + Pv o'Du)A{^{D_ r+„„ + _{D+ r_™) +A{R-P + P^~ 2M) {^{D_ r+„„ + r_™) 
- A - A _{D^ G,_ - _{D, A ^{D^ - „il?, A ^{D_ G„+} . (215) 

The wave operator on the left hand side of equation (|214p is (barring the —2A~^ term) A^'^q'D^"^ A^oVm A'"^ , which 
may be recast using the identity (with a an arbitrary constant) 

j^-2a ^ptm js^'^^D^ = A-'' oX>^'" QVra A"- - a{\ ^ R - P - 4M) + a^(\ - 2M) . (216) 

The identity (|216p . along with the operator equahties o^^^" o^^d = o^^^" o^'m and o^^^^ o2^+ = o^^^" o'^'- valid for zero 
spin, equations ([75)1 and (p5|) . brings the wave equation (|214p to a form similar to that already encountered for spin-±2 
and axial spin-0 modes, equations (|172p and (|199p : 

^-1/2 / ptS ^ ^pt± _ I _ + |p + |m) A-^l^k = . (217) 



8 4-"''4-''4^^ 

As written, the expression (j215p for the polar spin-0 source term S^^^ is valid for all harmonics, including dipole 

harmonics. For quadrupole and higher harmonics, I > 2, the tetrad connections in the source term S^^"^ can be replaced 
by spin-±2 components of the Weyl and Einstein tensors, in accordance with equations (|170p . valid in spherical gauge: 

2A 

r+™ + ^il?± r-™ = - - — Tvy— -r {±iV^ ±22?T C±2 + ^iV± ^2V± Cl^) , (218a) 

-uu +UU (( — lj(( + 2j 

+iV_ r+™ + _{D^ r_„„ = {i-i){i + 2) (+'^- +'^- ^++ + -'^+ ■ ^^^^^^ 

Physically, the right hand sides of equations (|218p are proportional to the polar parts of C±2 and G++ angularly 
lowered and raised to zero spin. 



44 



As discussed above, taken as a whole the wave equation (|217p is equivalent to a (coordinate and tetrad) gauge- 
invariant equation, being the difference of two gauge-invariant equations (a general equation and a monopole 
equation). However, the split between the left and right hand sides of equation (|217p is in general not gauge- 
invariant. In the case of the Schwarzschild geometry, where the background energy-momentum vanishes, and all 

(p) 

perturbations to the Einstein tensor are coordinate and tetrad gauge-invariant, the polar spin-0 source term Sq 
with the replacements ()218|) becomes a coordinate and tetrad gauge-invariant object. Therefore, at least in the case 
of the Schwarzschild geometry, and for quadrupole and higher harmonics, the left and right hand sides of the wave 
equation (I217P are individually equivalent to coordinate and tetrad gauge-invariant objects. 

The polar spin-0 wave equation (|217p has two unusual features compared to the other wave equations derived in 
this paper. The first is that (for quadrupole and higher harmonics) the spin-±2 Weyl tensor C±2 and its angular 

conjugate (7*2 7 in the combination given by equation ()218ap . appear in the source term Sq^\ equation ()215|) . That is, 

1 

the wave equation for polar spin-0 perturbation J\f appears to be sourced by spin-±2 modes of the Weyl tensor (even 

in the case of the Schwarzschild geometry) . 

The second unusual feature of the wave equation (|217p is that its homogeneous solutions decay with conformal 
1 1 
factor as A/" cx A~^''^ (so that the dimensional perturbation A~'^J^ with the same dimension as the spin-0 Weyl 

tensor Co goes as A~'^J\f (x A~'^/'^), which is a factor of A~^^^ faster than expected for the non-propagating spin-0 

companion of a propagating spin-±2 mode [such as the axial component Cq"', whose homogeneous solutions fall off 
as A~^, equation (|202p ]. This anomalous behavior is related to the presence of the spin-±2 Weyl tensor as a source 

of the wave equation. Since at least one of the spin-±2 components C±2 is propagating, it remains an active source 

1 

for the polar spin-0 perturbation Af even in a hypothetically flat empty region far from the black hole. Dimensional 

analysis indicates that the part of sourced by the driving spin-±2 Weyl component is expected to go with conformal 

1 , , _ _ 1 

factor A as A/" cx ASq cx A'^C±2, which according to the homogeneous solutions (|179p for C±2 goes as A/" ~ A^^ in 

1 

asymptotically flat empty space far from the black hole. Thus the sourced part of J\f is expected to dominate the 
homogeneous part far from the black hole. 

The homogeneous solutions of the wave equation (|217p are, similarly to equation ([75]) . 

Af = A-3/2e— *n^,(r) oY^e, 0) , (219) 
where the radial eigenmodes n^i{r) satisfy 

{-^-^' + [^(' + 1) - ^0 ■ (220) 

As discussed in the previous paragraph, in a hypothetical asymptotically flat empty region far from the black hole, 

1 

the behavior of the polar spin-0 perturbation Af is expected to be dominated not by a homogeneous solution (|219p , 

(p) 

but rather by the driving spin-±2 Weyl component in the source Sq . 

The polar potential [l^l 1) - | - |i? |P |M] A in the wave operator of the homogeneous polar spin-0 wave 
equation (|220p is not the same as the Zerilli polar potential even for the Schwarzschild geometry [g^, let alone 
for the Reissner-Nordstrom geometry [gl, [9^ . There is no contradiction because equation (|220p does not describe 
decoupled eigenmodes, but rather it describes the radial part of homogeneous solutions (I219P of a wave equation (I217P 
whose source sjf^ does not vanish even for the Schwarzschild geometry. 

The polar part of the spin-0 wave equation (|156ap in the Newman-Penrose formalism has (apparently) led to the 
wave equation (|220p , but it has not led all the way to an independent derivation of the Zerilli- Moncrief polar potential 
y*^"'"-' [gI, [9^ for decoupled gravito-elecromagnetic eigenmodes in either the Schwarzschild or Reissner-Nordstrom 
geometries. It would be good to accomplish such a derivation, but it may not be trivial to do so, and here I take the 
matter no further. 



VI. ELECTROMAGNETIC WAVES 



This section derives the equations that describe electromagnetic waves in perturbed self-similar black hole 
spacetimes. Subsections IVI Al and IVI Bl characterize perturbations of the electromagnetic field and the Maxwell's 
equations that govern their evolution. Subsequent subsections consider the propagating (spin-±l) components of 
electromagnetic waves fi jVI Cp . monopole electromagnetic modes ( WI Pp . and the axial f WI Ep and polar ( ^Vl Fp 
parts of the electromagnetic spin-0 wave equation. 
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A. Electromagnetic field 

The electromagnetic field i^„m is a bivector, an antisymmetric tensor. As earlier discussed in W Al with regard to 
the bivector structure of the Weyl tensor, this means that the electromagnetic field is a 6-component object with a 
natural complex structure. The real part of the eletcromagnetic bivector field is the electric 3-vector E, which changes 
sign under spatial inversion, while the imaginary part is the magnetic 3-vector B, which does not change sign under 
spatial inversion. 

The natural complex structure motivates writing Maxwell's equations as the single complex equation 

SkS? + ^e^;"") D'^Frnn = ^TTJi , (221) 

whose real (vector) and imaginary (trivector) parts represent respectively the source and source- free parts of Maxwell's 
equations. Absent magnetic monopoles, the current ji is purely electric, a real vector. Maxwell's equations (|22ip can 
immediately be recast as 

D'^Fki = 2Tr ji , (222) 

where the complexified electromagnetic field F is 

Fm = i (^C^r + y^r") (223) 

[as earlier in the definition (|13ip of the complexified Weyl tensor, the factor of 1/2 on the right hand side of 
equation (|223p is introduced because then the operator = \ {^T^i^ + t^fei™") is a projection operator, satisfying 
= P] In an orthonormal tetrad basis, the electromagnetic field Fki can be represented as a pair of electric and 
magnetic 3-vectors E and B 

F ^ {E B) . (224) 

The complexified electromagnetic field Fki then has the structure 

F = i ( 1 i){E + iB) . (225) 

Thus the magnetic part of the complexified electromagnetic field is i times the electric part. In tensor language, the 
complexified electromagnetic field satisfies the symmetry 

Fu = -£.,™"F,„„ , (226) 

as follows immediately from the definition (j223p . Thanks to this symmetry, the independent components of Fki 
constitute a single complex 3-vector, containing the electric field E and the magnetic field B in the complex 
combination E + iB. 

The complexified electromagnetic tensor Fki changes sign under spatial inversion, because spatial inversion not only 
transforms E —E and B B, but also changes the sign of the antisymmetric tensor Ekimm thereby effectively 
changing the sign oi i in E + iB . 

With respect to a general Newman-Penrose tetrad, the non-vanishing components of the complexified 
electromagnetic field Fki, and their expressions in terms of the electromagnetic field tensor Fki, a-nd the electric 
and magnetic fields E and B, are 

Fo = Fu. = \ (Fuv + F+-) = i (Er + iBr) , (227a) 

F±i = = F„+ = i [Eg + iBg T i{E^ + iB^)] . (227b) 

As with the Weyl tensor, the only non-vanishing components of the complexified electromagnetic field tensor Fki are 
those with bivector indices uv, h — , v+, or u-. Components Fki with bivector indices v~ or u+ vanish: 



(228) 
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In the unperturbed self-similar background, only the polar spin-0 component Fq of the electromagnetic field is 



non-zero, and it equals ^ the unperturbed radial electric field Er 

F„ = ^h = kk . (229) 

The unperturbed value of the electric field can be taken to define what is meant by the electric charge q interior to r 
in the unperturbed background, and the associated dimensionless charge Q = q/A: 



Fo = i q/A^ = \QIA. (230) 

Like the complexified Weyl tensor, equation (|138p . the complexified electromagnetic field tensor F^i can be resolved 
into polar (p) and axial (a) parts which respectively do not and do change sign when the azimuthal angular tetrad 
axis 70 is flipped in sign, equation ((54)) . 

Fs = /^i^^ + F'i"^ . (231) 

The angular conjugate of Fu is by definition 

F* = FjP) - F^f> . (232) 

The angular conjugates i^^; of the spin components of the complexified electromagnetic field tensor are, in terms of 
the electromagnetic field tensor F/j,;, and also in terms of the electric and magnetic fields E and B, 

n ^Wu.-F^_)^\ (Er - iBr) , (233a) 

FXi = [Ee - iBg ± i{E^ - iB^)] . (233b) 

Being tetrad- frame quantities, all components Fs of the complexified electromagnetic field tensor (and their angular 
conjugates) are coordinate gauge- invariant, although the split between the unperturbed and perturbed parts of the 
polar part of Fq, the only component whose unperturbed value is non- vanishing, is not coordinate gauge- invariant. The 
spin-0 component Fq (and its angular conjugate) are not only coordinate gauge-invariant but tetrad gauge-invariant 
with respect to all 6 arbitrary infinitesimal tetrad transformations, including those tetrad freedoms that are fixed in 
spherical gauge, equations (jl06p . The spin-±l components F^i (and their angular conjugates) are not tetrad gauge- 
invariant with respect to arbitrary infinitesimal tetrad transformations, being changed by a tetrad transformation 
that varies the antisymmetric part of Lpv+ (or its angular conjugate i/'d-). The following combinations of F^i or F\ 

and Fq are tetrad gauge- invariant, but not coordinate gauge-invariant: 

F±\ T {^v-v - ip+v)Fo , (234a) 

U— -U 

Fl^T hv- ~ip-v)Fo ■ (234b) 

The one component of the complexified electromagnetic field tensor that is coordinate and tetrad gauge-invariant 
and vanishing in the background is the axial spin-0 component Fq"'^ 

Fo"^ = 5(^0 - F*) = = ^iBr , (235) 

which is i« times the radial component Br of the magnetic field. 

The electromagnetic energy-momentum tensor is given by AnTmn = FmkFn^ — j'ymnFkiF''K The components of the 



electromagnetic energy-momentum tensor in a Newman- Penrose tetrad are given in the Appendix, equations (jA3p . 
To linear order, equations (|A3|1 reduce to 

AttTuv = 47rr+_ = = 2FoF* , (236a) 

4TTn+ = FyuFv^ = T {Q/A)F±i , (236b) 

U- UV u- 

AttTv- = FvuFv- =T(Q/A)Fri (236c) 
AttTvv = 47rr++ = . (236d) 

uu 

In the unperturbed background, the non-vanishing components of the electromagnetic energy-momentum tensor 
simplify to 

AnTuv = 47rT+_ = i(QM)2 . (237) 
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B. Maxwell's equations 



The central result of this subsection is the set of linearized Maxwell's equations (|241[) which describe the 
evolution of electromagnetic perturbations in spherically symmetric self-similar black hole spacetimes. The resulting 
electromagnetic wave equations (j248l) will be applied in subsequent subsections. 

In a general Newman-Penrose tetrad, Maxwell's equations (|222p are explicitly, 



- r+™ AF^i T o'^v Fq - iiD^p F^i = 271 A jv , 

-uu ^ u ^ u 

r+„+ AF^i T o2)±Fo-±i2)«F±i = 27rAj± , 

-u- V 

where the differential operators s3Dm are defined by 



uvv sr_+„ + (s + 2)r 



a± - s r_++ + s r„«+ + (s + 2) 



which reduce in the unperturbed background to 



(238a) 
(238b) 

(239a) 
(239b) 

(240) 



To linear order of perturbations on the self-similar background, Maxwell's equations (|238p reduce (in a general 
gauge) to 



T \ a'SvA-^Q T qVv Fo ~ F^i = 27rA> , 

^ u u u 

1 

In the unperturbed self-similar background, Maxwell's equations (|24ip simplify to 

T }^qVvA-^Q ^ 2TTAjv , 
T^o'DiA-'Q = . 



(241a) 
(241b) 



(242a) 
(242b) 



From Maxwell's equations (|222p and the antisymmetry of the electromagnetic field tensor Fm it follows that the 
electric current ji satisfies the charge conservation law 



D'ji = . 

In a general Newman-Penrose tetrad, the electric current conservation equation (I243P is explicitly 

A^D^ji = - +iK^jv - -iK A]u + -1®'+ Aj_ + Aj^ = , 
where the differential operators sSJ^ are 



2)' 



A^ 
A^ 



dv -I- s Ft, 



(s + 1) F_+„ + F_„+ + (s + 2) F+„_ 



9±-sF_++-f (s + 1) 



A-1 



which reduce in the unperturbed background to 



Linearized, the electric current conservation equation (j244p reduces to 

- A - Aju + -iV+Aj_ + ^{D_ A jV = 



(243) 
(244) 

(245a) 
(245b) 

(246) 
(247) 



The pair of Maxwell's equations (|241ap can be combined into wave equations for the complexified electromagnetic 
tensor Fs in much the same way as was done earlier for the complexified Weyl tensor Cs , equations (jl65p . Taking one of 
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the operators sS'„ of equations (|245p times equation (|241a[) minus another of the operators s23'„ times equation (j241bp . 
the choice of operators being guided by the electric current conservation equation (|247p . yields equations that look like 
sourced wave equations for each of the (perturbations of the) spin components Fg of the complexified electromagnetic 
tensor (in a general gauge): 

(oI?t%2?s + o2?^%2?±)^o = -i±iI)^(oD«A-iQ±4^Aj.) +i±iI?^(o3±^"'Q±47r^j±) , (248a) 
(^^iV^^^ ±iVu + ±iV^^ ±iV^)f±i = ^o'3L{±o^vA-^Q + AnAjv)-^oVv{±Q^±A-'Q + A7:Aj±y (248b) 

As with the wave equations (|165p for the Weyl tensor Cs , the wave equations (|248p for the electromagnetic tensor Fg 
are not independent wave equations: the spin components Fs are not independently adjustable, but rather fluctuate 

in harmony with each other in accordance with Maxwell's equations (|24ip . 

1 

Equation (|248ap appears to constitute two separate equations for Fq^ but the two equations are in fact equivalent. 
The mathematical derivation of the equivalence is almost identical to that which demonstrates the equivalence of the 

two wave equations (|165ap for Cq, except that it is necessary to use the electric current conservation equation (|247p 
instead of the Weyl current conservation equations (|164p . 

It should be remarked that the electric current conservation equation (|247p is insufficient to imply the definition (|245p 
of the operators sD'^ for arbitrary spin s, in particular for the zero spin operators o2)m that appear in the wave 
equation (j248bp for F±i. Actually, it is fine to replace the differential operators s'^'m on the right hand sides of 
both wave equations l|248p with their unperturbed limits ±iDm, since the combined quantities in parentheses that 
they operate on are already of first order (although the separate terms inside the parentheses are not individually of 
first order). The advantage of retaining iiS'^ in equation (|248aP is that it makes application of the electric current 

V 

conservation equation (|247p . as in the previous paragraph, transparent. As regards equation (j248bp . the advantage of 
qI)'^ as given by equations (I245p . besides the resulting symmetrical relation between the expressions (|162l) and (|245p 
for s^'m ^-nd s^'m the expressions (|154p and (|239p for and s^m, is that this choice climates terms equal 
to first order quantities times (d )A and (d± + 2T-t-vu)A that would otherwise appear on the last 

"7/, III' ^ lilt' 



-u+' 



line of equation (jA7p . and that would then need to be converted to (unperturbed) currents using equation (|242p . 
The operators o®ln defined by equations (|245p are, modulo the conformal factors A, the same as those used by 
Chandrasekhar before equation (204) on page 238. 

The tricky parts of the derivations of the wave equations for the spin-±l components F^i in WI C[ and for the 
axial spin-0 component F^^ in WI El are the parts that depend on the unperturbed dimensionless interior charge 
Q on the right hand side of equations (|248p . In the case of this tricky part of the right hand side of the wave 
equation (j248bp is given as equation (jA7p in the Appendix, which linearizes (in a general gauge) to 

{qT)'^qT>v-o'Z)'^q^±)a-^Q = 2q{2AC±i^A±2'DuA-^T+vv ^ ±2'D^T^v+) - AT qVu A'^ Q , (249) 

the last term on the right hand side of which can be converted to unperturbed electric currents j; using equation (|242p . 

In the case of the axial spin-0 component Fq"''' , the tricky part of the right hand side of the wave equation (|248ap 
is given as equation (jA8p in the Appendix, which linearizes (in a general gauge) to 

oS.-+i2)'_ oS+-+iS)L* o2):++iS'* oK) ^^'Q - 4 q a c^^'^^-a [(r_,+-r+,_) oi?„+(r+n--r_„+) oA,] a-'q , 

(250) 



the last set of terms on the right hand side of which can be converted to unperturbed electric currents j; using 
equation (I242p . 



C. Spin-one (propagating) components of electromagnetic waves 

The spin-±l components F±i of the complexified electromagnetic field are of particular interest because they 
describe the propagating components of electromagnetic waves. The main result of this subsection is the sourced 
spin-±l electromagnetic wave equation (|251l) . whose homogeneous solutions (|255p have radial part satisfying the 
generalized spin-±l Teukolsky equation (|256p . 
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If the unperturbed black hole is uncharged, Q = 0, then the spin-±l components F±i are both coordinate and tetrad 
gauge-invariant and vanishing in the background, and therefore represent real physical perturbations. However, in 
the more general case where the black hole is charged, Q ^ 0, then the spin-±l components Fj-i are not tetrad gauge- 
invariant, being changed by an infinitesimal tetrad transformation of the antisymmetric part of ipv+- For Q ^ 0, it is 

in principle possible to use the tetrad gauge freedom to set F±i — 0. Indeed, in application to the Reissner-Nordstrom 
geometry, Chandrasekhar [l^ page 240 adopts precisely such a gauge, calling it the "phantom gauge" . However, 
from a physical perspective it is unnatural to adopt a gauge where the propagating component of a field is set to 
zero. As noted previously, expression (|234ap . the combination F±i =p (ipv+ — (p+v)Fq is tetrad gauge-invariant. In 

situations where the propagating components of the electromagnetic field are much larger than the static component, 
S> 1^0 1 J which an observer would interpret confidently as light waves, setting F±i to zero would require setting 
the supposedly small Lorentz transformation ipv+ — f+v to some large value, which would be quite unnatural. Thus 

the spin-±l components F^i of the electromagnetic field have physical content notwithstanding the tetrad gauge 
ambiguity, and it makes sense to consider the wave equation that governs their propagation. 

The linearized wave equation for the spin-±l components F^i of the electromagnetic field is, from equations (I248b|) 
and ((249)) . 



(iipt- ^^Vu + ±il?^^ ±iI?tJF±i = 5±i , (251) 
where the source term 5±i is, in a general gauge, 

5±i ^±Aq(2C^i~ + A'^T -f 27ryl2 r+™ ju + 27rfo£i± A> - oI?« A j±) . (252) 

V V ^uu -u-J ~uu V \ ^ u u ) 

In spherical gauge, where r+u^ vanishes, the source term equation (|252p . reduces to 

— u — 

5±i =±AQf2C'±i-±2l?«^~^r+.™) -H27r^2p °^_^27rfoD'±^> -o2?«v4j±) . (253) 

V V -mm/ -mm V \ ^ u u J 

In spherical gauge, and for non-dipolc modes, I 7^ 2, the connection F+mu is determined by the spin-±2 components 

— MM 

C±2 of the complexified Weyl tensor, equation (|168ap . Equation (|253p for 5±i thus shows that electromagnetic waves 

are sourced by two separate effects, first by gravitational waves, the first sets of terms on the right hand side of 



equation (j253p (including the term proportional to the unperturbed electric current jj), and second by perturbations 

to the electric currents j; , the last set of terms on the right hand side of equation (|253p . If the unperturbed black 



hole is uncharged, Q = [which among other things implies that the unperturbed electric current vanishes, ji — 0, in 
accordance with equations (|242p ]. then the gravitational wave source vanishes, but if the black hole is charged, then 
gravitational waves are an inevitable source of electromagnetic waves. 

Much of the discussion in W CI concerning the propagating components C±2 of gravitational waves, starting from 
the paragraph containing equation p78p and on to the end of the subsection, carries through with little change for the 
case of the propagating components F±i of electromagnetic waves. The angular conjugates F'^-y satisfy the angular 
conjugate of equation (|25ip 



{Du + TiI?T± ^^2?^ = . (254) 



The angular Hermitian operators +12?^ and _il'^'''_i2?+ inside the wave operators of equation (I25ip and its 

angular conjugate (|254p have the same eigenvalues (and complex conjugate angular eigenfunctions), and it follows 
that the temporal-radial (i.e. non-angular) part of the eigenfunctions of the wave operators are the same for both F^\ 
and its angular conjugate F^^. 

The wave operator on the left hand side of equation (|25ip is identical to the wave operator discussed in ijlllEl 
equation ([75]) . for the case of spin one, s = 1. The homogeneous solutions (when the source terms vanish) of the wave 
equation (I25ip and its angular conjugate (|254p are, similarly to equation ([78)) . 

|f } = (r/r)^^/^A-2e--,i/Kr) I (255) 

where ±ifi{r) depends only on conformal radius r, and satisfies the second order ordinary differential equation 

{- ^ + 0^ ± 2k)' + [1(1 + 1) - 1 - 6C] a| = , (256) 
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which agrees with the generahzed Teukolsky equation ([T]) given in the Introduction, since, equation 

-6C = -2(i?- P) - 2Pl + 6M . (257) 

Equation (j256p has two hnearly independent eigensolutions, which may be identified as ingoing and outgoing based 
on their properties in a hypothetical asymptotically flat region at large radius, as already discussed in ijlllFl In a 
hypothetical aysmptotically flat region at large radius, the radial eigensolutions sfuji (with s = ±1) go over to the 
asymptotic solutions sZi discussed in §111 F[ equation (|94|) . Because the wave operators for modes of opposite spin 
are complex conjugates of each other, equation (j256p . it follows that the ingoing and outgoing radial eigensolutions 
sfuji of opposite spin can be taken to be complex conjugates of each other, g/^" = (_s/^/^*) , in complete analogy to 
equation ([5^ . The propagating components of Fg (x A~'^e~'^'^* sfuji are those that fall off most slowly at large radius, 
which identifies the positive spin component F^i as describing the propagating component of ingoing electromagnetic 
waves, and the negative spin component F_i as describing the propagating component of outgoing electromagnetic 
waves: 

F^\ cx A~'^e'''^\ifJi - ^-2^*g-it^(t+r-*) (propagating, ingoing) , (258a) 
cx A-^e-"^'.if°f ~ J^-2^*^-^u.{t-r') (propagating, outgoing) . (258b) 

The other of the two modes for each spin corresponds to the non-propagating, short-range spin-(— s) partner of a 
propagating mode of spin s: 

cx A-'^e-'"^^_Jl^i - A-2(r*)-ie"*"(*+''*' (non-propagating, ingoing) , (259a) 
F°f cx A-^e''"^\J°Y' ~ A-^{r*y^e-'"^^''''"> (non-propagating, outgoing) . (259b) 

D. Electromagnetic monopole modes 

For monopole modes, / = 0, all spin-±l quantities vanish identically, and therefore the spin-±l components F±i of 
the complexified electromagnetic field tensor cannot describe monopole modes. The bottom line of this subsection is 
the unsurprising result that the equation (j266p for the polar part of the electromagnetic monopole mode simply enforces 
conservation of interior electric charge, and does not describe any electromagnetic wave. The axial electromagnetic 
monopole mode vanishes identically, equation (|268p (at least in the absence of magnetic monopoles). 

An equation for the monopole follows from the linearized Maxwell equation (j241a[) . which can be written (in a 
general gauge, and for arbitrary harmonics) 

A'^(dv+T^y++r+y^A-^Q + 2AoVvF^'''^ ± 2A±iI?^ F±i = AQ(r_„+ - r+«_) :f 4tt A^ jv , (260) 

where Q is the generalized interior electric charge defined by (notation: script Q is the generalized interior electric 
charge, to be distinguished from the italic Q dimensionless unperturbed interior electric charge related to the 
unperturbed electric charge q hy Q = q/A) 

Q = 2A2F(|p) , (261) 

whose unperturbed value is just the electric charge q interior to radius r in the unperturbed background, 
equation (j230p . 

Q^q . (262) 

The generalized interior charge Q is tetrad gauge-invariant, but it is not coordinate gauge-invariant, and moreover 
the split between its unperturbed and perturbed parts is not coordinate gauge-invariant. However, the unperturbed 



part Q = q, which transforms under a coordinate gauge transformation as 

Q^Q + A{e^d, + e"9„)Q = Q " 4^^(e"L - e"j„) , (263) 
does become coordinate (and tetrad) gauge-invariant in any region where the unperturbed background electric charge 



and current vanish, — 0, such as in a hypothetical asymptotically flat empty region far from the black hole. In such 
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a region where the unperturbed background charge and current vanish, the unperturbed Maxwell's equations (|242|) 



imply that the unperturbed interior electric charge Q goes over to a constant, the charge g, of the black hole, 

n 

Q ^ q, = constant , (264) 

a coordinate and tetrad gauge-invariant quantity. 

For the monopole mode, I = 0, all non-zero spin objects in equation (|260p vanish identically [including T_y_^ — T_^y_ 
and r+u_ — r_„+, whose expressions ()104c|) in terms of vierbein components ipmn contain only spin-±l terms], and 
equation (|260p simplifies to 

A^(dv+T^v++T+v-)A-^Q + 2AQVvF^''^ ^T'iTrA^jv (1 = 0). (265) 

The polar part of the electromagnetic monopole equation p65p is, in a general gauge, 

A^(dv+T-v++T+y-)A-~^Q = T'i'^A^jv (1 = 0). (266) 

In spherical gauge, where equation ()112ap holds, equation (I266P simplifies to 

a„g = =F47rA2> (/ = 0). (267) 

u u 

Equation (I266|) . or its spherical gauge version (|267p . simply expresses conservation of (the monopole part of the) 
generalized interior electric charge Q to linear order. 

The axial part of the monopole equation (|265p is, in a general gauge, qDv Fq"'' = 0, which integrates to 

Ffj") = (1 = 0) , (268) 

that is, the monopole axial spin-0 component of the electromagnetic field tensor is identically zero, a coordinate and 
tetrad gauge-invariant statement. 



E. Axial spin-zero component of electromagnetic waves 

The previous two subsections, WI C I and I VI D I have given a complete description of electromagnetic perturbations: 
WI CI covered spin-±l electromagnetic waves, which encompass all electromagnetic perturbations except monopole 
modes; while §VI CI treated the special case of monopole modes. Thus in principle there is no need to probe the 
electromagnetic equations any further. 

However, for the reasons enunciated at the beginning of W Gl it is desirable to consider wave equations for the 
polar and axial parts of the spin-0 component Fq of the electromagnetic field tensor, and since the Newman-Penrose 
formalism does provide a spin-0 wave equation for electromagnetic perturbations, equation (|248ap . it seems a good idea 
to explore where this equation leads. This subsection gives the wave equation (|269p for axial spin-0 electromagnetic 
modes, while the next subsection WIF I considers the wave equation for polar spin-0 electromagnetic modes. 

The axial spin-0 component Fg°'\ equation (|235p . of the complexified electromagnetic field tensor, which equals 
times the radial component Br of the magnetic field, is the only component of the electromagnetic field tensor 
that is coordinate and tetrad gauge-invariant and vanishing in the background (an assertion that includes the 

propagating spin-±l components, which are not tetrad gauge-invariant). The linearized wave equation for i^g"^ 
is, from equations (|248ap and (|250|) . 

(opt" ^Vr, + oI?t± Ffl'') = 5^°) , (269) 
where the source term 5q°'' is, in a general gauge, 

5^"^ = -Q^c^'^^-i-7rA2[(r_„+-r+„_)i„-(r+„_-r_„+)j,] -7rA(_iP+jL-+iP_jV) . (270) 

Axial gauge, equation (|197p . simplifies the the source term Sq""^ for axial spin-0 electromagnetic waves, as was 

previously found in the source term S'q"^ for axial spin-0 gravitational waves. In axial gauge, the source term iSq"'' 
reduces to 

S^,"^ =~QACt^ -7tA{_iV+j_-+,V_u) . (271) 
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Equation (|27ip for iSq"-* shows that the axial spin-0 component of electromagnetic waves is, like the spin-±l 
components, sourced by two separate effects, first by gravitational waves, the first term on the right hand side of 
equation (|27ip , and second by perturbations to the electric currents ji , the last set of terms on the right hand side of 
equation ((?7T|) . 

Much of the discussion in W Gl concerning the axial spin-0 component Cq"-* of gravitational waves, starting from 
the paragraph containing equation (j202|) and on to the end of the subsection, carries through with little change for 

the case of the axial spin-0 component Fq"''' of electromagnetic waves. The homogeneous solutions (when the source 



terms vanish) of the wave equation ()269|) are, similarly to equation ([78|) 

i^o^"^ = A-2e— *o/i:^ (r) oil™ {e, cf>) , (272) 

where the radial function afj^"^ (r) in equation (j272p depends only on conformal radius r, and satisfies the second order 
ordinary differential equation 



ofLf = . (273) 



The scalar potential l{l + 1)A in equation (|273p agrees with that given by [i^ [his equations (91)-(93) for s — 1]. 

Equation (|273p has two linearly independent eigensolutions, which may be taken to be complex conjugates of each 
other, and which may be identified as ingoing and outgoing based on their properties in a hypothetical asymptotically 
flat region at large radius, as discussed in mil Fl In a hypothetical aysmptotically flat region at large radius, the 



radial eigensolutions o/^l; go over to the asymptotic solutions qZi discussed in ^III F[ equation ([87ll . 

°feout ° ' (^r ) ^ e as r* ^ oo . (274) 



F. Polar spin-zero component of electromagnetic waves 

The main result of this subsection is the sourced wave equation ()279p for the polar spin-0 component of 
electromagnetic waves. The polar wave equation (|279p is the same as its axial counterpart (|235p . aside from the 
difference in source terms. 

The polar spin-0 component Fq^^ — ^{Fq + Fq) — ^E^ of the electromagnetic field, which equals i the radial 
component Er of the electric field, is both coordinate and tetrad gauge- invariant. If the unperturbed black hole 
is uncharged, Q = 0, then Fq is also vanishing in the background, and its perturbation then represents a real 
physical perturbation. However, in the more general case where the black hole is charged, Q ^ 0, the polar spin- 
electromagnetic field Fq^^ does not vanish in the background, and then the split between its unperturbed and 
perturbed parts is not coordinate gauge-invariant. Indeed, as long as the unperturbed spin-0 electromagnetic field 

1 

-Fq^' is non-constant, it is possible in principle to eliminate the perturbation Fq^' altogether by a coordinate gauge 
transformation arranged to satisfy 

- A{e^d, + e^du)^^ = , (275) 

1 

in much the same way that it is possible in principle to eliminate the perturbation Cq^' to the polar spin-0 Weyl tensor, 
equation (|209p . However, as with the Weyl tensor, there is reason to expect that the polar part of the electromagnetic 
spin-0 wave equation (j248ap is not vacuous, but genuinely encodes the evolution of polar spin-0 electromagnetic 
perturbations, as is certainly the case if the unperturbed black hole is uncharged, Q = 0. 

In a sense it is straightforward to obtain a wave equation for the polar spin-0 electromagnetic perturbation: the 
Maxwell equation (|260p is (coordinate and tetrad) gauge- invariant (that is, the split between the left and right 
hand sides of the equation is not gauge- invariant, but the equation taken as a whole is gauge- invariant), and the 
monopole equation (j265p is similarly gauge-invariant, and therefore the difference between the general and monopole 
equations is gauge-invariant, and moreover vanishing in the background. When this equation is combined with 
equation (|276p below, the result will be the desired electromagnetic polar spin-0 wave equation, which as a whole 
will be gauge-invariant (though the split between left and right hand sides will in general not be gauge-invariant). 
However, subtracting the monopole equation does not eliminate the contribution of the monopole component of the 
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electromagnetic field from the wave equation, because the wave equation involves perturbed derivatives dmQ of the 
unperturbed monopole, which in general do not vanish because in general the unperturbed monopole is not constant. 

In the case of gravitational polar spin-0 waves, i 3VH| it was found most natural to cast the wave equation in 
terms of the polar spin-0 quantity N defined by equation (|190p . which had the property that its unperturbed part 
was a constant, equation l|19ip . a coordinate and tetrad gauge-invariant object. For electromagnetic polar spin-0 
waves, there does not appear to be a polar spin-0 quantity whose unperturbed part is always constant. The closest 
analogous quantity is the generalized interior electric charge Q, equation (j26ip . whose unperturbed part is constant, 

_____ 

equation (|264p . over any region where the unperturbed background electric charge and current vanish, j; — 0, such 
as in a hypothetical asymptotically flat empty region far from the black hole. 

The first Maxwell equation (|241ap was already recast in terms of the generalized interior electric charge Q in 
equation (j260p . Similarly recasting the second Maxwell equation (I241bp in terms of Q gives (in a general gauge, and 
for arbitrary harmonics) 

d±Q + 2A oPi^i"^ ± 2 A ±iI?„F±i = 2Q d^A - 2AQ r+™ t 471^^ j± . (276) 

Combining the polar parts of equations (j260p and (j276p yields the following prototype wave equation for the generalized 
interior charge Q (in a general gauge) 



+1^'^ A{d, + r_,+ + r+„_) + A{du + r^u- + r_„+) _ a-\iV_ 9+ - d 



Q 

= 2A-^Q{^{D_ a+ A + a_ A) - Q{^iV_ r+,„ + r_„„ + +iV_ r+„„ + r_„„) 

- A7r{^,^'^Aj,~ _,^'^A,U) , (277) 

where the differential operators s'^'m ^-re given by equations (I245|) . Equation ()277p has not yet had its unperturbed part 
subtracted. In spherical gauge, where equations (|112ap and (|212p hold, and r+t,„ — T+uv according to equation (|109ep . 

-uv -vu 

the prototype wave equation (|277p becomes 



= 4A-1QAA - 2Q{^iV_ r+™ + _i2?+ r_„„) - 4^(+iD:, A j„ - Aju) , (278) 



which again has not yet had its unperturbed part subtracted. It is seen that if the unperturbed part Q of the 



generalized interior electric charge were constant, so that all its derivatives vanished, dmQ ~ 0, then Q on the left 

1 

hand side of equation (|278p could be replaced by its perturbation Q. However, the unperturbed interior electric 



charge Q is in general not constant. Therefore, when the unperturbed part of equation (|278p is subtracted, the 



residual contributions from perturbed derivatives of Q must be taken over to the right hand side, where they become 

incorporated into the source term. 

1 

Acting on the perturbed part Q of the generalized interior electric charge, the wave operator on the left hand side 
of equation (|278p is o'D^"^ o'Dm A~^ . This, along with the operator equalities o'D^^ qT^v — o^^^" o^'u and qI^^'^ oT^+ = 
o'D''~ o'D_ valid for zero spin, brings the electromagnetic polar spin-0 wave equation (|278p in spherical gauge to the 
form 

(opt" + oPti >(f ) = 5^^) , (279) 



in which the perturbation Q of the generalized interior electric charge has been replaced by the equivalent perturbation 

■Ap) - 1 

~ 2" 



-^cT^ ~ 5^ of the polar spin-0 electromagnetic field, so as to bring out the similarity between the polar spin-0 



wave equation (|279p and its axial counterpart (|269p . The polar spin-0 source term Sq^'' on the right hand side of the 
wave equation (|279p is, in spherical gauge, 

5*^^ = A-'QM- r+,„ + _iV+ r_„„) - TT{+iVuAj^ - AjJ 

- \{+iVuA-^k + _iV,A-^du-A-\iV_d+-A-^_iV+d^)Q, (280) 
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10 

where the perturbations dm Q to the directed derivatives of the unperturbed interior electric charge Q are 

1 

d,n Q = {'f mu dv + if mv 9„) Q , (281) 



which would vanish if the unperturbed interior electric charge Q were constant. 

The expression (|280|) for the source term iSp^^ in spherical gauge is valid for all harmonics, including monopole 

and dipole harmonics. For quadrupole and higher harmonics, I > 2, the combination ^iD_T_^yu + -I'D+^-vu of 
(p) 

tetrad connections in Sq can be replaced by the combination of spin-±2 Einstein components G++ and G__ given 

by equation (j218bp . valid in spherical gauge. Equation (|280p for S^f^ shows that, like the other components (spin-±l 
and axial spin-0) of electromagnetic waves, the polar spin-0 component of electromagnetic waves is sourced both by 
gravitational waves, the first two sets of terms on the right hand side of equation (|280p . and by perturbations to 
electric currents, the last two sets of terms on the right hand side of equation (j280p including the terms involving 



derivatives of the unperturbed interior electric charge Q. As with the other spin components of electromagnetic waves, 
if the unperturbed black hole is uncharged, (5 = 0, then the gravitational wave contribution to the electromagnetic 
source vanishes, but if the unperturbed black hole is charged, Q ^ 0, then gravitational waves are an inevitable source 
of electromagnetic waves. 

The electromagnetic polar spin-0 wave equation (I279|) is the same as its axial spin-0 counterpart (|269p apart from 

the difference in the polar and axial source terms S(^^ and Sq°'\ Therefore the homogeneous solutions for the polar 

1 

spin-0 perturbation Fg^-* are essentially the same as those (|272p for the axial spin-0 component Fq"'\ The polar 

radial eigenfunction aflj'i^i'r) satisfies the same differential equation (|273p as its axial counterpart, and its asymptotic 
behavior is likewise the same (|274p as its axial counterpart. 

VII. SCALAR WAVES 

This section derives the equations that describe massless scalar waves in perturbed self-similar black hole spacetimes. 

A. Massless scalar field 

No massless scalar field is known in nature. Nevertheless, massless scalar fields are of interest because they describe 
hypothetical waves of spin zero that propagate at the speed of light, and as such provide a theoretical counterpart 
to the known waves of spin one (electromagnetic) and two (gravitational) that propagate at the speed of light. The 
wave equation for a massless scalar field ■0 is the massless Klein-Gordon equation 

D'^dm^ = . (282) 

Without loss of generality, the scalar field 4' can be taken to be real in an orthonormal tetrad. In quantum field 
theory, it is natural to consider a complex scalar field ■0, in which the field V' and its complex conjugate ip* describe 
particles and anti-particles respectively, but this is not needed for the present paper. 

The scalar field is by assumption a scalar, a coordinate and tetrad gauge-invariant object. However, the split 
between the unperturbed and perturbed parts of the scalar field is not coordinate gauge-invariant. 

Being scalar, the field is unchanged by angular conjugation, which flips the sign of the azimuthal tetrad axis 7^, 
equation l|54p . The scalar field tp is thus a purely polar quantity. For a scalar field there only polar waves, no axial 
waves. 

The energy-momentum tensor of the scalar field is given by AnTmn = T^m'^n — ^ImnT^kT^^ , page 483 of (60| . where 
7r„i is defined to be the gradient, or momentum, of tp, 

-Km = dm'lp ■ (283) 

The components of the scalar energy-momentum tensor in a general Newman-Penrose tetrad are given in the Appendix, 
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equations (jA4p . To linear order, equations (|A4p reduce to 



47rT+_ = TTyTiu , (284a) 
/^t:Tvv = TT^TTt, , (284b) 

ii-jf U U 



4TrTy+ = nvTT± , (284c) 



47rT„_ = ^i^TT^ , (284d) 
47rr„„ = 47rT++ = . (284e) 



B. Massless scalar waves 

The main result of this subsection is the sourced massless scalar wave equation (|292p . whose homogeneous 
solutions (|293p have radial part satisfying the generalized spin-0 Teukolsky equation (|294p . 
In a Newman-Penrose tetrad, the massless scalar wave equation (|282p is explicitly 

{^i^Ad, + -iS; Aa„ - +iS)'_ Ad+~ = , (285) 

where the differential operators are defined by equation (|245p . The wave equation (|285p linearizes to 

{+iT)'^Ad, + Adu - +iV_ Ad+ - _iV+ Ad_)^ = . (286) 
If the unperturbed wave equation is subtracted from the linearized wave equation (|286p . and the residual contribution 



from perturbed derivatives of the unperturbed field taken over to the right hand side, then the result is a wave 

1 

equation for the perturbation ijj (it is worth remarking here that d'^'^dm = A^^ g^?^™ A^ o^-'m A^"^) 

opt™ o-Dm A-^ = 2a , (287) 

with source a [the factor of 2 on the right hand side of equation (|287p is inserted so that it disappears from 
equation (I292p below] given by the following combination of derivatives of the unperturbed scalar field: 

<J = ^iVuAd^- _iV^Adu + +iV_Ad+ + _iV^Ad_ 

+ A^ [- {du - Luu + r+«- + r»«+)<9« - {k - r„™ + r_.+ + r+.„)a„} , (288) 

1 

where the perturbations dm "0 to the directed derivatives of the unperturbed scalar field are 

10 

5™ V = {'fmu dy + iprav 9„)V' ■ (289) 



If the unperturbed scalar field -0 is constant, as for example occurs asymptotically in a hypothetical asymptotically 
flat empty spacetime far from the black hole, then the source cr, equation (|288p . vanishes identically, and the wave 
equation (|285p is (coordinate and tetrad) gauge-invariant. However, in the more general case that the unperturbed 



scalar field ^ is non-constant, then the split between the unperturbed and perturbed parts of the scalar field is not 
coordinate gauge- invariant, and the wave equation (j287p is then likewise not coordinate gauge-invariant. Indeed, 

as long as the unperturbed scalar field is non-constant, it is in principle possible to eliminate the perturbation ?/; 
altogether, by a coordinate gauge transformation arranged to satisfy 

Af^dJ^ = . (290) 

1 

However, if the perturbation i/; to the scalar field varies on scales much smaller than the scale over which the 

unperturbed field "0 is varying, then the gauge choice (j290p can be accomplished only by making the supposedly 
small coordinate perturbation e™ large. Thus the gauge choice (|290p is not natural, and it is increasingly less natural 
for perturbations of higher frequency in regions where the unperturbed scalar field is more slowly varying. 
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It can be argued that a more natural gauge choice is simply to demand that the source term a, equation (|288p . 
vanishes 

a = 0, (291) 

a condition that is automatically satisfied in any region where the unperturbed scalar field is constant, such as in 
a hypothetical asymptotically flat empty region far from the black hole. However, in the general case that the 
unperturbed scalar field is not constant, the gauge condition (|29ip represents a rather complicated condition on the 
vierbein perturbations ipmm which may not be ideal for numerical calculations. Moreover, in the asymptotically flat 
region far from the black hole, the question of how the condition (|29ip translates into a condition on the vierbein 
perturbations ipmn may become numerically delicate, because the gauge condition (|29ip is satisfied automatically 
rather than being a gauge condition. 

On balance, it may be preferable to adopt whatever gauge happens to be most numerically convenient, and to 
accept that the scalar source term (|288l) may be complicated. 

The identity (|216p . along with the operator equalities o'D^'" ^Dy = qD^^ qV^ and o2^^^ 0^+ = o'D^^ (yD_ valid for zero 
spin, equations (|75p and brings the linearized scalar wave equation (|287p to a form similar to those previously 
encountered for gravitational waves, equations (|172p . (|199p . and (|217p . and electromagnetic waves, equations (|25ip . 
and (17751) . 

" ti n 1 

oPt" qVu + oV'<^ aV^ + \{-R + P + 2M) A-^ilj = a . (292) 

The wave operator on the left hand side of equation (I292p differs from the wave operator discussed in §111 El only 
by the addition of the terms ^(— i? + P + 2M), which are functions only of conformal radius v. The homogeneous 
solutions of the scalar wave equation (|292p (when the source term vanishes) are, similarly to equation ([78]) . 



^ = A-^e-'^''^^i{T)Yi^{B, 0) , (293) 
whose radial part ipujiir) depends only on conformal radius r, and satisfies the second order differential equation 

^--u;^+ [1(1 + 1) - R + P + 2M]A^Mr) =0 , (294) 

which agrees with the generalized Teukolsky equation ([1]) given in the Introduction. The scalar potential [l{l + 1) — 
R + P + 2M ]A i n equation ((294)) agrees with that given by ^ [his equations (91)-(93) for s = 0]. 

Equation (l294p has two linearly independent eigensolutions, which may be taken to be complex conjugates of each 
other, and which may be identified as ingoing and outgoing based on their properties in a hypothetical asymptotically 
flat region at large radius, as discussed in mil Fl In a hypothetical aysmptotically flat region at large radius, the 
radial eigensolutions ip^^i go over to the asymptotic solutions qZi discussed in i jlllFl equation (187)) . 



out _ .out/, , „^c.r•• as r* oo . (295) 



VIII. SUMMARY 



This paper has presented the foundations of the theory of perturbations of spherically symmetric self-similar black 
holes. The theory encompasses the limiting cases of the Schwarzschild and Reissner-Nordstrom geometries, and of 
more general stationary spherical spacetimes. No additional condition on the nature of the energy-momentum has 
been imposed beyond those of spherical symmetry and self-similarity (for example, the energy-momentum is not 
required to be an ideal fluid). 

The primary motivation has been to explore the mass inflation instability inside black holes, where the geometry 
departs hugely from any stationary geometry, but given the length and technicality of the present paper, it has seemed 
best to defer application, including all issues of numerical computation, to subsequent work. 

Wave equations have been derived for the propagating components of gravitational, electromagnetic, and scalar 
waves, and it has been shown how the non-propagating, short-range components are related to the propagating 
components. The special cases of the (non-propagating) monopole and dipole components of gravitational 
perturbations, and monopole components of electromagnetic perturbations, have been treated. The homogeneous 
solutions of the wave equations are separable with respect to a suitable set of conformal coordinates: a conformal 
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|. From the outset, each vierbein perturbation is 



time t and conformal radius r, and the usual polar and azimuthal angular coordinates and (p. A generalization of 
the Teukolsky equation is given as equation ([1]) in the Introduction. 

Although separable, the wave equations are not decoupled. In the presence of non-vanishing background energy- 
momentum (that is, except in the limiting case of the Schwarzschild geometry), waves of each spin - tensor (spin-2), 
vector (spin-1), or scalar (spin-0) - are sourced by perturbations of other kinds, so that waves of different spin are 
inevitably coupled to each other, and convert from one kind to another. This is unlike the case of the Fricdmann- 
Robertson- Walker metric of cosmology, where spatial translation and rotation symmetries ensure that tensor, vector, 
and scalar modes are decoupled completely from each other P, 0j_ ( 

The paper follows the Newman-Penrose null tetrad formalism [33, 
required to be expandable in spin- weighted harmonics of definite spin (as opposed to the more general system of tensor 
harmonics), and it is shown that this requirement can be imposed consistently without any gauge conditions. If each 
vierbein perturbation is expanded in spin-weighted harmonics, then it automatically follows that each component 
of the tetrad-frame Riemann tensor, and consequently Einstein and Weyl tensors, is expandable in spin-weighted 
harmonics of definite spin. The spin of each perturbation follows the general rule that it equals the sum of the -l-'s 
and — 's of its covariant angular indices in the Newman-Penrose formalism. 

Although the Newman-Penrose wave equations resolve into spin- weighted harmonics with no gauge conditions, if in 
addition it is required that each tetrad connection (Ricci rotation coefficient) be expandable in spin- weighted harmonics 
of definite spin, then this imposes a set of 8 gauge conditions, here termed "spherical" gauge, equations (|106p . Spherical 
gauge proves to provide a particularly nice gauge for expressing the source terms of propagating gravitational waves, 
and also for exploring polar (but not axial, apparently) spin-0 perturbations, including monopole modes. 

The Newman-Penrose formalism provides wave equations not only for propagating components of waves, but also 
for non-propagating components. Wave equations for polar and axial spin-0 modes are derived in this paper, partly in 
an attempt to make some contact with the traditional description of per turbations to the Schwarzschild and Reissner- 
Nordstrom geometries in terms of polar and axial perturbations [isl. Il9l. [6^. [s^. [o^ [osj . The Newman-Penrose spin-0 
wave equation does correctly reproduce the Regge- Wheeler axial potential for the Schwarzschild geometry, but this is 
the only potential that it recovers without further manipulation. The Newman-Penrose formalism yields not decoupled 
wave equations for polar and axial spin-0 perturbations, but rather sourced wave equations. It would be useful in 
future work to see how the Newman-Penrose spin-0 wave equations for case of the Reissner-Nordstrom geometry 
might be decoupled so as to recover the Zerilli-Moncrief polar and axial potentials. 
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APPENDIX A: NEWMAN-PENROSE COMPONENTS 

This appendix gives expressions for the Newman-Penrose components of various quantities, in the notation of the 
present paper. 
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Riemann tensor 



The Newman-Penrose components of the Riemann tensor Rkimn hi terms of the tetrad connections Tkmn and their 
directed derivatives are, from equation (|37p . 
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2. Weyl currents 



The Newman-Penrose components of the Weyl currents Jimn on the right hand sides of equations (|153p for the 
evolution of the complexified Weyl tensor are 
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3. Electromagnetic energy-momentum tensor 



The Newman-Penrose components of the energy-momentum tensor Tmn — (l/47r)(i^„ifeF„'^ — \lmnFkiF^'') of the 
electromagnetic field are given by: 
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4. Scalar energy-momentum tensor 

The Newman- Penrose components of the energy-momentum tensor Tmn — (l/47r) (7rm7r„ — ^7mn7rfe7r'^) of a massless 
scalar field "0, with gradient (momentum) 7r„j = dmip, are given by: 
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5. Some operator commutators 

The following equation is referred to in WBl equation (|166p : 
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The following equation is referred to in WIB| equation (|249|) [compare equation (jASp . which is quite similar]: 
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The following equation is referred to in WIB[ equation (|250p [compare equation (|A6p . which is quite similar]: 
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